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1.  PRODUCTION 

Nucleation  and  condensation  phenomena  are  of  fundamental  importance 
in  many  fields  such  as  atmospheric  physics,  astrophysics,  aeronautical 
sciences,  nuclear- reactor  technology  and  the  material  sciences.  They  must 
be  considered  in  weather  prediction,  aerosol  formation,  metallurgical  tech- 
niques, high-speed  wind  tunnel  design,  steam-turbine  nozzle  design,  rocket 
performance  and  high- intensity  molecular  beams.  The  processes  include 
condensation  of  supersaturated  vapour  with  or  without  foreign  nuclei 
crystal  formation  from  vapour,  melted  and  solid-liquid  solutions,  phase 
changes  in  solid  states,  crack  nucleation  in  solids,  and  crystallization 
of  polymers . Consequently , nucleation  and  condensation  problems  have 
been  investigated  intensively  in  recent  years  in  diverse  fields.  The 
results  of  these  studies  were  summarized  by  Hirth  and  Pound  ( 1Q631 
Zettelmoyer  (1969,  1977)  and  Abraham  (1974).  ’ 

The  formation  and  behaviour  of  condensed  matter  are  controlled  by 

occur>  in  “any  cases,  in  heterogeneous  phases. 
The  micrqphysical  mechanism  by  which  nucleous  atoms  or  molecules  become 
activated  and  grow  to  larger  size  clusters  is  still  poorly  understood, 
especially  for  heterogeneous  nucleation  and  condensation.  Although  homo- 
geneous nucleation  and  condensation  are  not  frequently  found  in  nature  and 
not  very  applicable  to  practical  problems,  the  concepts  established  through 
the  study  of  homogeneous  nucleation  and  condensation  provide  a basis  for 
understanding  the  physically  more  complicated  process  of  heterogeneous 
nucieation  and  condensation.  Much  effort  has  been  expanded  by  many  resear- 

hal^a^^century1^*1^^1^  understandlng  of  homogeneous  nucleation  for  over 

Although  it  was  an  ambitious  attenpt  to  construct  a theory  for 
a nonequilibrium  process  like  nucleation  and  condensation,  the  early 
investigators  only  managed  to  produce  an  approximate  theory.  Nevertheless, 
it  succeeded  in  obtaining  remarkable  agreement  with  available  experimental 
data.  The  theory  connects  equilibrium- rate  distribution  of  clusters  of 
near- critical  size  to  the  steady-state  distribution  of  clusters.  A very 
small  cluster  is  assumed  as  a spherical  drop  of  liquid.  The  free  energy 
of  the  cluster  is  assumed  to  be  equal  to  the  free  energy  of  the  number  of 
molecules  in  the  bulk  liquid  plus  the  free  energy  due  to  the  surface  of 
the  cluster.  The  supersaturated  vapour  is  viewed  as  a mixture  of  ideal  gas 
species,  each  species  corresponding  to  a subcritical  cluster  of  a given 
size.  The  equilibrium  distribution  of  these  clusters  is  then  calculated. 

th  an  equilibrium- rate  theory,  the  rate  of  nucleation  can  be  predicted  to 
be  proportional  to  the  equilibrium  concentration  of  the  cluster  of  critical 
size  at  which  the  free  energy  of  the  cluster  is  a minimum. 

The  crudeness  of  the  calculation  of  the  equilibrium  distribution 
was  recognized  and  improved  using  the  more  detailed  methods  of  statistical 
mechanics.  However,  the  classical  theory  still  serves  adequately  in  the 
interpretation  of  available  experimental  results  except  for  the  case  of 
very  high  cooling  rates  as  found  in  supersonic  free-molecular  jets.  Since 
classical  theory  is  based  on  a quasi- steady  state  nucleation  concept,  it 
becomes  inapplicable  when  the  time  required  to  achieve  steady- state 
nucleation  is  large  compared  with  a characteristic  time  governing  the  change 
of  state.  In  ordinary  cloud  chambers  and  nozzle  flows,  the  characteristic 
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time  of  nucleation  is  generally  too  short  to  be  accounted  for.  In  rapid 
expansions  of  free-molecular  and  cluster  beams  it  takes  a longer  time  for 
the  vapour  to  adjust  to  a steady- state  than  to  pass  through  a nozzle. 
Further,  classical  theory  assumes  the  free  energy  of  formation  of  a 
critical  droplet  can  be  expressed  in  terms  of  the  macroscopic  surface 
tension  of  the  bulk  liquid  and  its  volume  free  energy.  In  highly  super- 
saturated situations,  the  critical  nucleus  may  be  exceedingly  small, 
containing  few  molecules.  Consequently,  the  macroscopic  concepts  of 
continuum  thermodynamics  become  invalid. 


In  order  to  under stand  the  nucleation  mechanism  of  small  clusters, 
theoretical  studies  were  made  from  the  standpoint  of  irreversible  statisti- 
cal mechanics  by  assuming  intermolecular  interactions  between  atoms  or 
molecules  and  using  molecular  dynamics  or  Monte-Carlo  methods.  In  order 
to  solve  the  nonequilibrium  process  of  nucleation,  the  kinetics  of  cluster 
formation  were  examined  theoretically.  In  addition,  the  early  stages  of 
clusterization,  especially  dimerization,  were  investigated  experimentally. 

Nucleation  and  condensation  processes  in  ccmpzaBsible  flows  are 
of  importance  in  the  aerospace  sciences  in  designing  supersonic  and  hyper- 
sonic wind  tunnels  and  rocket  nozzles.  In  turn,  such  flows  provide  funda- 
mental data  needed  for  an  understanding  of  nucleation  and  condensation 
processes.  The  heat  release  due  to  condensation  can  be  a significant  flow 
disturbance.  Consequently,  the  equations  of  motion  with  heat  addition  are 
solved  and  compared  with  experimental  data  to  test  the  validity  of  different 
nucleation  theories.  In  this  respect,  well-established  flows  in  cloud 
chambers,  supersonic  nozzles  and  shock  tubes  have  provided  useful  data  in 
understanding  macroscopic  processes  of  nucleation  and  condensation. 

Reviews  of  such  investigations  were  made  by  Wegener  (1969,  1975  (a),  (b) , 
1977) • 
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Experiments  with  cluster  beams  in  free-molecular  jets  have  provided 
means  of  studying  the  material  properties  of  small  clusters.  They  also 
yielded  fundamental  insight  into  the  kinetic  process  of  nucleation.  Most 
experiments  on  gas-phase  nucleation  in  cloud  chambers,  continuum  nozzle 
flows  or  rarefaction- wave  flows  were  limited  to  finding  the  super saturation 
ratios  at  which  condensation  occurs  during  an  expansion.  A quantitative 
interpretation  of  such  experiments  is  hampered  by  the  difficulty  in 
accounting  for  small  complex  processes  in  terms  of  a simple  result,  say,  in 
an  integral  form.  To  provide  all  the  desired  information  about  condensation, 
the  ideal  experiment  would  be  the  one  in  which  the  stepwise  event  from 
monomer  to  droplet  could  be  observed.  Molecular  beam  sampling  of  a free 
jet  is  capable  of  providing  such  information.  A number  of  such  studies 
have  been  done  in  this  decade. 


Fogs  and  clouds  are  very  important  phenomena  of  nucleation  and 
condensation  from  a meteorological  point  of  view.  Similar  phenomena  on 
a smaller  scale  are  experienced  daily  closer  to  home  by  using  sprays, 
kettles,  cigarettes  or  pipes.  However,  the  mechanisms  of  nucleation  and 
condensation  in  these  cases  are  substantially  different  from  those  found 
in  supersonic-nozzle  flows,  expansion  waves  and  free- Jet  flows.  Nucleation 
and  condensation  occur  more  heterogeneously.  Liquid- soluble  substrates 
or  insoluble  liquid- wet  substrates  mostly  dominate  the  nucleation  process. 
Physical  and  chemical  surface  characteristics  of  the  substrates  such  as 


hydraphi id. city,  hydrqphobicity,  contamination,  and  surface  catalyticity  are 
much  more  important . Since  the  droplet  size  is  usually  of  the  order  of 
1 ~ 100  jam,  condensation  is  affected  by  more  macroscopic  features  of  the 
system  such  as  turbulent  coalescence,  turbulent  miirlng  and  precipitation. 
Droplet  growth  is  affected  predominantly  by  the  factors  of  continuum  flow 
rather  than  those  on  the  molecular  level.  Their  features  are  too  different 
to  discuss  them  simultaneously. 

In  the  present  review,  the  subject  matter  is  restricted  to  nonequi- 
librium flows  with  nueleation  and  condensation.  Special  emphasis  is  given 
to  the  theories  of  homogeneous  nueleation  and  heterogeneous  nueleation.  The 
cited  references  are  limited  to  those  published  since  1965.  Unfortunately, 
it  was  hot  possible  to  include  a comprehensive  review  of  the  Russian  litera- 
ture. Perhaps  this  will  be  done  by  Russian  researchers  in  this  field. 

It  should  be  noted  that  no  claim  is  made  for  cocprehensiveneas . 

The  field  is  just  too  large  for  that.  Our  modest  aim  was  to  provide  a 
starting  point  for  new  researchers  in  this  area.  Neither  is  any  claim  mad** 
for  a complete  critical  analysis  of  deficiences  in  the  various  nueleation 
and  condensation  theories.  Certainly  the  use  of  a correction  factor  for 
surface  tension  or,  condensation  coefficient  ac  and  others  such  as  r,  to 
try  to  fit  ejqserimental  data  is  a very  unsatisfactory  procedure.  Surer 
foundations  will  have  to  be  laid  in  the  future. 


2.  NUCLEA1I0N 

A vapour  consists  of  clusters  of  molecules  bound  together  by  their 
intermolecular  interactions.  Nueleation  from  the  vapour  phase  can  be 
considered  as  a kinetic  process  that  vapour  molecules  and  their  clusters 
combine  with  each  other  to  form  larger  clusters  (dimers,  trimers, 
i-mers) . The  elementary  process  may  be  expressed  as  a series  of  stepwise 
reactions  consisting  of  bimole eular  combination  and  unimole cular  dissocia- 
tion. 


A1  +A1^A2 
*2  + A1  A3 


A.  , + 


Ai  + 


where  A^_  means  a single  vapour  molecule  (monomer)  and  A^  is  a cluster 
containing  i molecules  of  Ai  (i-mer) . Clusters  Ai  (i  = 2,  3,  ...)  grow 
into  Ai+i  by  addition  of  a single  molecule  through  bimole cular  combination 
(condensation)  and  decay  into  Ai_i  by  loss  of  a single  molecule  (evaporation) 
(Fig.  l).  In  addition  to  these  processes,  there  may  exist  other  processes 
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such  as  the  formation  of  large  clusters  by  the  collision  of  two  small 
clusters  (Ai  + Aj  -» Ai+j) , the  dissociation  of  large  clusters  into  smaller 
ones  (Ai+j  -+  Ai  + Aj) , or  three  body  reactions  (Ai  + Aj  + A*  -+Am  + An) . 
However,  these  processes  can  be  excluded  as  being  too  infrequent  except 
the  initial  dimerization.  It  is  recognized  that  the  initial  dimerization 
involves  a third-body  reaction  (see  Chapter  4.4)  such  as 

Al+\*Al^iA2*Al 

^ ^ + B Ag  + B 

where  B means  a single  molecule  of  the  third  body  (e.g.,  carrier  gas) . The 
phase  change  of  condensation  takes  place  when  the  growth  process  of  combine- 
tion  dominates  the  decay  process  of  cluster  dissociation,  thereby  yielding 
increasing  numbers  of  large  clusters.  Thus,  if  this  multistate  kinetics 
of  clusterization  is  known,  it  is  possible  to  calculate  the  time  history 
of  the  concentration  of  clusters  according  to  the  rate  equations  subject 
to  the  initial  conditions  of  clusters  and  the  physical  conditions  of  the 
system. 


When  only  bimolecular  reactions  are  taken  into  account,  the  number 
density  of  clusters  of  size  i,  nj.,  is  given  by  solving  the  following  equations: 


dn. 

d r = 1i~  Ii+1  (i  ^ 2) 


d*1! 

dt“ 


(2) 


where  Ii  is  the  rate  of  formation  of  the  i-mer  given  by 


^i  °i-lni-l  ” eini 


(3) 


where  Ci_i  is  the  rate  at  which  molecules  condense  on  a cluster  A^_^  and 
ei  is  the  rate  at  which  molecules  evaporate  from  a cluster  Ai  (Fig.  1). 

The  vapour  molecules  strike  a cluster  Ai  of  'i  collision  cross-section 
Si  with  the  mean  molecular  speed  Ur,  and  only  a certain  fraction  ac  of  / 
these  molecules  sticks  to  the  cluster.  The  condensation  rate  ci  is 
expressed  as 

ci  ■ Wl  • "l 

According  to  the  kinetic  theory  of  molecular  collision,  the  mean  molecular 
speed  can  be  given  by  (8kT^-»  i)l/2  where  T is  the  temperature , k the 
Boltzmann  constant,  and  mi  i tne  reduced  mass  of  the  cluster  Ai  and  monomer 
Ai  (=  mi/(i+l),  m the  mass’of  monomer).  The  collision  cross-section  Si  is 
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given  by  vrj2  where  r*  is  the  mean  collision  radius  of  and  Ai.  Thus, 
when  the  vapour  is  an  ideal  gas  (ni  * Py/kT),  the  condensation  rate  is 
rewritten  as 


°1  * “o  4*ri‘ 


}1kr 


Two  particular  states  of  the  system  can  he  specified;  the  steady 
case  in  which  the  concentration  r.j.  does  not  change  with  time,  dn^/dt  = 0; 
and  the  quasi-equilibrium  case  in  which  the  rate  of  cluster  formation  is 
zero,  Ii  = 0,  although,  strictly  speaking,  an  equilibrium  state  does  not 
exist  in  a nucleation  process.  In  the  latter  case,  from  Eq.  (3) 

e.  n . _ 

i _ i-l»e 


where  the  subscript  e denotes  the  equilibrium  state.  Small  deviations  of 
the  ci_i  and  ei  from  the  equilibrium  values  are  assumed,  c^_]_  e « ci_i, 
ei,e  "*  ®i*  Thus,  to  predict  the  reaction  rates  ci  and  ei,  it’is  required 
to  know  the  statistical  thermodynamic  properties  and  the  reaction  cross- 
sections  as  a function  of  the  relative  kinetic  energy  of  the  reactants, 
their  mutual  orientation  and  the  internal  energy  of  rotation  and  vibration. 
Even  with  such  knowledge  available,  it  would  be  too  difficult  to  derive  an 
analytical  formulation  of  a general  nature  for  the  exact  calculation  of  the 
nucleation  process  except  for  dimer  or  trimer  formation.  For  the  latter 
cases,  with  steady  progress  in  the  understanding  of  intermolecular  forces, 
it  appears  possible  to  dfevelop  such  a rigorous  kinetic  theory.  It  Is  for 
these  reasons  that  almost  all  existing  (practical)  nucleation  theories 
bypass  considering  intermediate  microscopic  stages  and  adopt  a thermo- 
dynamic or  statistical  approach  in  which  only  the  initial  and  final  stages 
are  quoted  without  time  dependence. 


2.1  Statistical 


c Theo 


The  statistical  thermodynamic  theory  of  nucleation  is  a steady- 
state  theory  in  which  the  concentrations  of  clusters  are  assumed  to  remain 
in  a steady  state: 


(i  > 2) 


with  which  Eq.  (2)  gives 


I.  = const  ■ I 


(i  >2) 


To  maintain  the  system  strictly  in  such  a steady  state,  clusters  must  be 
removed  from  the  system  when  they  grow  up  to  a certain  terminal  size  i^, 
ni*  = 0,  and  monomers  are  added  to  the  system  at  the  rate  i„  I.  From  Eqs.  (3) 
and  (5),  the  rate  of  cluster  formation  I is  related  with  the  steady-state 
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I 


concentration  of  clusters  n^  and  the  equilibrium  concentration  of  clusters, 
ni,e> 


I = c.n.  (i  .2>i> 

1 i,e  V“i>e  l»  - - 


1) 


'i-1  "i,e  n. 


1-1, e 


(i=ij 

Sunning  the  above  equations  divided  by  cj.ni  e leads  to 


(7) 


Here,  the  reciprocal  of  the  steady- state  rate  of  cluster  formation  is  equal 
to  the  sum  of  the  reciprocals  of  all  the  equilibrium  condensation  rates  of 
clusterization.  This  is  a fundamental  equation  of  nucleation  rate  which 
has  been  discussed  by  many  investigators  (Zettelmoyer,  1969,  1977).  For 
clusters  of  large  size,  the  discrete  variables  i in  Eq.  (7)  are  replaced 
by  differentials  to  convert  the  summation  to  an  integral  form. 


Evaluation  of  the  nucleation  rate  from  Eqs.  (7)  or  (8)  requires 
knowledge  of  the  condensation  rate,  , and  the  equilibrium  concentration 
of  clusters,  nijC.  The  former  earn  be  predicted  by  Eq.  (4)  from  the  mole- 
cular kinetics.  By  introducing  the  free  energy  of  formation  of  clusters 
Af,  the  concentration  of  clusters  in  thermodynamic  equilibrium  can  be 
expressed  in  the  form 


n 


i,e 


(9) 


The  free  energy  of  cluster  formation  is  equal  to  the  elevation  of  the  free 
energy  for  a transition  from  monomers  to  the  cluster, 


i 


(10) 


where  Gi  and  Gi  are  the  free  energy  of  monomer  and  cluster,  respectively. 
The  free  energy  of  the  cluster  can  be  determined  a priori  by  confining 
quantum-mechanical,  statistical -mechanical,  and  Monte-Carlo  techniques, 
or  it  can  be  approximated  in  terms  of  the  properties  of  the  bulk  material, 
being  a function  of  cluster  size. 

Usually,  the  free  energy  of  cluster  formation,  ZGi>  increases 
with  the  size  of  cluster,  and  is  a maximum  kt  some  value  of  i,  say  i*, 


6 


at  which 


(11) 


The  cluster  corresponding  to  the  maximum  value  ZG^  is  called  "the  critical 
nucleus  or  critical  embryo",  and  has  the  lowest  equilibrium  concentration 
(Fig.  2).  Once  clusters  achieve  this  critical  size,  the  addition  of  another 
molecule  to  the  cluster  causes  a decrease  in  the  free  energy  of  cluster 
formation.  This  means  a high  probability  that  they  will  grow  to  macroscopic 
size  with  relatively  unhindered  conditions.  The  value  of  £G^  does  not 
change  very  much  with  respect  to  i around  i = i* . Evaluating  the  integral 
in  Eq.  (8)  by  expanding  ZGi  about  its  maximum  with  respect  to  i leads  to 
an  expression  of  I*  in  the  form 


I*  ■ I 


* 

c.  n, 
i l,e 


(12) 


where  Zc  is  known  as  the  Zeldovich  factor  and  given  by 


(13) 


The  nucleation  rate  given  by  Eq.  (7)  also  can  be  reduced 
as  Eq.  (12)  where  the  Zeldovich  factor  is  given  by 


Z 

c 


n 


i,e 


to  the  same  form 


(13a) 


The  expression  of  nucleation  rate  by  Eq.  (12)  means  that  the  steady- state 
rate  of  nucleation  is  equal  to  the  equilibrium  value  with  a correction 
factor  Zc  which  takes  the  value  of  10“1  ~ 10" 2 . 


The  equilibrium  concentration  in  Eq.  (8)  is  dominated  by  the 
free  energy  of  cluster  formation,  ZGi*,  which  depends  on  the  statistical 
thermodynamic  properties  of  clusters.  Thus,  the  nucleation  rate  is 
predominantly  controlled  by  this  free  energy.  The  quantitative  estimation 
of  the  free  energy  of  cluster  formation  is  the  major  problem  in  quasi- 
steady nucleation  theory.  The  physical  reasoning  and  mathematical  simpli- 
fications involved,  such  as  nonequilibrium  factors,  choice  of  cluster 
partition  function,  re- evaporation,  etc.,  provide  a variety  of  additional 
effects.  Recently,  such  statistical  thermodynamics  of  nucleation  was 
generalized  by  Derjaguin  (1974) . 


2.2  Homogeneous  Nucleation 

The  simplest  description  of  the  free  energy  of  cluster  formation 
is  the  liquid-drop  model,  often  called  the  capillarity  approximation,  which 
has  been  used  successfully  for  many  years  in  the  various  areas  of  nucleation 
and  condensation.  The  liquid  drop  model  is  based  on  a simple  physical 
picture  of  the  cluster.  It  is  assumed  that  the  cluster  is  viewed  as  a 
uniform  spherical  drop  of  liquid,  which  is  basically  bulk- like  except  for 
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the  presence  of  the  bulk- like  surface.  Thus,  the  free  energy  of  formation 
of  a droplet  composed,  of  i monomers  is  assumed  to  be  equal  to  the  sum  of 

(i)  the  change  in  the  free  energy  of  i molecules  from  the  vapour  into 
the  condensed  phase,  and 

(ii)  the  free  energy  of  the  generating  surface  on  the  bulk  liquid  in  an 
amount  equal  to  the  surface  area  of  the  droplet, 


**i  = “ Hv)  + Sjo  (l4) 

where  p./  and  py  are  chemical  potentials  per  molecule  in  the  equilibrium 
bulk- liquid  phase  and  the  supersaturated  vapour  at  the  pressure  prevailing 
outside  of  the  drop,  respectively;  a the  macroscopic  interfacial  free 
energy  existing  between  the  bulk  liquid  and  the  vapour  (the  surface  tension) , 
surface  area  of  a sphere  having  a volume  equal  to  i times  the  molecular 
volume  Vjh£  in  the  bulk  liquid  phase, 

Si  . 4^  = 4y  ( fa  1 Vml  ) 2/3 

where  r is  the  radius  of  the  drop.  A more  generalized  theory  for  non- 
spherical  and  anisotropic  clusters  was  studied  by  Ziabicki  (1968a,  b) . 

The  critical  cluster  is  given  by  Eq.  (ll)  with  Eq.  (l4)  as 


r*  = 


By  - Bj 


I*  = 5l r*2 

3 By  - Bi 


which  gives 


iC«  « ^ 0T»2  * ( ■ 

3 3 \BV-B, 


mi 


(15) 


(16) 


where  the  superscript  * means  the  critical  nucleus.  The  difference  in 
chemical  potentials  is  expressed  as 


= kT  in  s (17) 


where  s is  the  super  saturation  or  the  ratio  of  the  pressure  of  the  conden- 
sable vapour  py  to  the  vapour  pressure  of  the  condensed  phase  (saturation 
pressure)  ps,  or  s = Pv/ps.  From  Eq.  (12)  with  Eqs.  (4)  and  (16) , the 
nucleation  rate  is 


<18> 


8 


This  is  the  nucleation  rate  of  the  so-called  classical  ( Becker-Doring) 
nudeation  theory.  The  condensation  coefficient  is  usually  hut  unjusti- 
fiably assumed  as  ac  = 1. 

The  Zeldavich  factor  is  given  by  Eqs.  (13)  and  (l4) 


The  factor  usually  is  of  the  order  of  10“^  ~ 10"2,  being  relatively  insig- 
nificant  compared  with  the  remaining  pre-exponential  value  in  Eq.  (18)  of 
about  1025 . a more  exact  prediction  for  this  factor  was  studied  by  White 
(1969)  using  a steady-state  random  walk  analysis.  The  usual  expression  for 
Zc,  Eq.  (19) , was  found  to  be  quite  accurate. 

The  applicability  and  validity  of  the  classical  theory  of  nuclea- 
tion are  confined  by  the  assumptions  of  steady-state,  bulk  thermodynamic- 
properties  and  large  cluster  size.  In  ordinary  cloud  chambers,  supersonic 
nozzle  and  shock  tube  flows,  the  relaxation  time  required  to  establish  the 
steady- state  concentration  distribution  of  clusters  is  of  the  order  of 
microseconds,  while  the  characteristic  time  of  change  of  the  flow  state 
extends  to  milliseconds.  Here  the  steady-state  assumption  should  be  valid. 
On  the  other  hand,  in  supersonic  free-molecular  jets,  to  be  discussed  in 
Section  4.4,  the  characteristic  flow  time  is  extremely  small  compared  with 
the  build- up  time  of  steady- state  nucleation.  For  such  a situation,  nuclea- 
tion must  be  treated  as  a completely  nonequilibrium  process  by  solving  a 
set  of  reaction  equations  (see  Sections  2.4  and  4.4). 

In  Eq.  (18),  the  concentration  of  monomers  participating  in  the 
nucleation  and  condensation  is  assumed  to  be  constant.  However,  in  the 
process,  the  vapour  molecules  are  consumed  owing  to  nucleation  and  conden- 
sation and  in  time  distributed  to  the  associated  clusters,  or  diffused 
away  from  the  region  of  interest  under  nonuniform  conditions  in  space. 

These  effects  were  studied  by  Katz  et  al  (1966) , Katz  and  Blander  (1973) > 
Heist  et  al  (1976),  Ziabieki  (1977),  and  Becker  and  Reiss  (1976,  1977). 

In  the  associated  vapour  of  distributed  clusters,  larger  supersaturations 
are  required  for  the  same  nucleation  rate  than  that  of  the  vapour  containing 
monomers  only,  owing  to  the  thermodynamic  properties  of  the  vapour  rather 
than  the  kinetic  formation  process.  Under  the  usual  conditions  in  diffusion 
cloud  chambers,  the  diffusion  effect  due  to  nonuniformities  does  not  cause 
serious  error. 

In  the  exponential  of  Eq.  (18),  (in  s)2  appears  through  r*2. 

Thus,  the  nucleation  rate  varies  quite  rapidly  with  supersaturation . It 
is  this  rapid  variation  of  the  nucleation  rate  which  gives  the  "critical 
supersaturation" , below  which  one  can  see  nothing  occurring  and  just  above 
which  nucleation  suddenly  occurs.  The  nucleation  rate  also  depends  very 
strongly  on  the  surface  tension  as  q3j  a 10$  change  in  a leads  to  changes 
in  the  nucleation  rate  by  ten  orders  of  magnitude.  Thus,  the  correct 
values  of  these  bulk  properties  are  very  important.  Mast  of  them  are 
readily  obtainable  for  common  vapours  from  handbooks  or  with  suitable 
extrapolations.  The  major  uncertainty  stems  from  the  possible  effect  of 
drop  size  on  the  value  of  surface  tension  when  the  size  is  small.  The 


« 


dependence  of  the  surface  tension  on  drop  radius  is  expressed  by  functioned, 
forms  like 


O’ 


where  Ca,  is  the  surface  tension  acting  on  the  flat  surface  of  the  condensate. 
However,  there  is  no  rigorous  analysis  which  is  suitable  for  condensation 
nuclei  of  100  or  fewer  molecules  of  the  complex  liquid.  Usually,  the 
macroscopic  flat-film  surface  tension  is  employed  for  nucleation  rate 
calculation.  An  improved  expression  for  the  nucleation  rate  with  macro- 
scopic surface  tension  as  a function  of  droplet  radius  was  developed  by 
Derjaguin  and  Prokhorov  (1974) . 

Fran  statistical  mechanics  the  concen- 
chemical  equilibrium  is  given  by 

<21> 

where  Z±  is  the  partition  function  for  a cluster  of  size  i and  p is  the 
chemical  potential  of  the  monomer  in  the  vapour . When  the  motion  of  the 
cluster  is  separated  into  translation  of  the  whole  cluster,  rotation  of 
the  whole  cluster,  and  internal  motion  of  the  molecules  within  the  cluster, 
the  partition  function  can  be  written  as 

Zi  ■IZtr«n8(1>Zrot(1>J>Zint<1^) 
i 

where  j represents  distinguishable  geometrical  configurations  of  the  cluster. 
The  translational  partition  function,  Ztrans,  is  independent  of  the  cluster 
conf iguration . The  rotational  partition  function,  Zr0£,  depends  on  the 
configuration  only  through  the  asymmetricity  of  moments  of  inertia.  The 
internal  partition  function,  Zint>  requires  complete  knowledge  of  the 
cluster  configuration.  When  Eq.  (21)  is  compared  with  Eq.  (9) > it  is  seen 
that  the  free  energy  of  cluster  formation  can  be  expressed  in  a different 
form  from  that  of  the  drop  model,  Eq.  (l4) . There  are  two  approaches  to 
improving  the  drop  model  in  this  way.  The  first  is  a hybrid  approach  which 
assumes  basically  the  drop  model  with  statistical  mechanical  corrections. 

The  second  is  to  evaluate  directly  the  statistical  mechanical  formulation 
for  the  free  energy,  calculating  the  partition  functions  of  clusters  for 
all  sizes  which  requires  both  a knowledge  of  the  interaction  potential  of 
molecules  of  the  cluster  and  suitable  techniques  of  calculation  such  as 
molecular  dynamics  and  Monte  Carlo  methods. 

The  first  approach  was  initiated  by  G.  M.  Pound  and  his  coworkers 
(the  so-called  Lo the- Pound  theory)  and  by  H.  Reiss  and  his  coworkers  (the 
so-called  Reiss-Katz-Cohen  theory)  (Zettlemoyer,  1977).  The  clusters  in 
the  vapour  may  behave  dynamically  like  large  gas  molecules,  having  free 
translational  and  rotational  motions.  These  motions  contribute  to  the 
partition  function  of  the  cluster  by  factors  which  are  equal  to  the  trans- 
lational and  rotatiocLpartition  functions  as  an  almost  rigid  body,  q.trans 


Statist! cal- thermodynamic  theory: 
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and  qpotj  respectively.  In  order  that  the  cluster  has  such  motions  in 
the  vapour,  certain  motions  of  the  cluster  in  the  liquid  hulk  phase  must 
he  deactivated  as  it  leaves  the  liquid  bulk,  since  the  total  nuniber  of 
degrees  of  freedom  of  the  cluster  is  conserved.  The  six  degrees  of  freedom 
used  in  qtrans  8:1(1  9rot  can  1)6  compensated  by  introducing  a "replacement 
factor"  qrep  which  is  related  to  the  internal  partition  functions  of  the 
cluster  (Lothe  and  Pound,  1968) . This  leads  to  an  increase  in  the  pre- 
exponential of  the  nucleation  rate  given  by  the  drop  model  in  the  form 


I = 


<4. 


rans 


qrot 


rep 


Lc  | 


(22) 


where  Ic^  is  the  nucleation  rate  given  by  Eq.  (18).  Assuming  that  the  six 
degrees  of  freedom  of  the  cluster  and  the  average  degrees  of  freedom  of 
the  bulk  liquid  have  the  same  partition  function,  Lothe  and  Pound  (1966) 
estimated  qtrans  'IrotArep  **  101?  for  water  vapour  nucleation.  In  view  of 
the  agreement  with  experimental  result,  this  inconsistency  was  called  the 
"translation-rotation  paradox"  by  Reiss  and  Katz  (1967) . They  considered 
that  the  centre  of  mass  of  a stationary  drop  fluctuates  with  high  probability 
over  a volume  comparable  to  its  physical  volume . A partition  function  for 
this  replacement  motion  nearly  cancels  the  partition  function  for  the  motion 
of  the  centre  of  mass  of  the  cluster  in  the  vapour.  Assuming  the  free 
energy  of  a drop  created  by  removing  a volume  with  fixed  boundaries  from 
the  bulk  liquid,  Reiss,  Katz  and  Cohen  (1968)  and  Lin  (1968)  showed  an 
increase  in  the  nucleation  rate  of  103  ~ 10°  for  water,  which  decreases 
the  predicted  supersaturation  only  by  8 ~ 10$.  Extensive  discussions  have 
appeared  on  this  matter  by  Lothe  and  Pound  (1966) , Reiss  and  Katz  (1967), 
Reiss,  Katz  and  Cohen  (1968),  Lin  (1968),  Parlange  (1968) , Stillinger  (1968), 
Abraham  and  Pound  (1968),  Lothe  and  Pound  (1968) , Cohen  (1970) , Reiss  (1970, 
1977),  Kikuchi  (1969,  1971,  1977),  Nishioka  et  al  (1971),  Nishioka,  Lothe 
and  Pound  (1972),  Blander  and  Katz  (1972)  and  Abraham  (1974a),  Reiss  (1977), 
Nishioka  and  Pound  (1977),  Nishioka  (1977),  Miyazaki  et  al  (1977).  However, 
because  of  the  current  disagreement  for  a physical  model  of  a cluster  there 
is  an  extreme  range  in  the  value  of  qtrans  (lrot/(lrep  frcan  10°  to  10^7 . 


To  resolve  the  disagreement  between  the  experiments  and  the  Lothe- 
Pound  theory,  Abraham  (1969,  1974a)  discussed  the  liquids  studied  in  the 
experiments  and  features  of  their  bulk  surface  structures  such  as  liquid 
bonding  characteristics,  surface  entropy,  molecular  dipole  moment,  acentric 
factors,  and  catalytic  effects  for  condensation  in  ions.  Hydrogen-bonded 
liquids  are  categorized  into  two  classes  of  oriented  polar  molecules  and 
non- oriented  quasi- spherical  molecules  at  the  liquid  surface.  The  former 
liquids  follow  the  classical  nucleation  theory,  whereas  the  latter  liquids 
follow  the  Lothe-Pound  theory.  He  suggested  that  the  failure  of  the  Lothe- 
Pound  theory  to  explain  the  nucleation  of  polar  liquids  was  due  to  its 
failure  to  account  for  the  deviation  of  critical  nucleus  properties  from 
bulk- liquid  surface  properties.  If  the  dipole- dipole  interaction  energy 
of  the  highly  oriented  surface  molecules  is  corrected  for  surface  curvature, 
the  surface  free  energy  of  the  drop  is  appreciably  greater  than  that  pre- 
dicted from  the  plane  surface. 

The  classical  model  including  the  modified  one  mentioned  above 
has  many  merits.  It  avoids  the  formidable  problem  of  treating  each  cluster 
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on  a molecular  basis  and  relates  the  free  energy  of  formation  to  measurable 
bulk  properties.  However,  these  properties,  such  as  surface  tension  and  surface 
area,  have  a questionable  interpretation  for  clusters  with  less  than  100  mole- 
cules and  in  addition  the  classical  model  has  concept! onal  difficulties  in  its 
inability  to  uniquely  relate  terms  in  the  statistical  model  to  terms  in  the 
classical  model  (Eqs.  9 and  21).  These  difficulties  can  be  avoided  by  a direct 
treatment  of  statistical  mechanics  of  small  clusters.  This  is  the  second 
approach  to  improving  the  drop  model.  For  this  purpose,  a static  technique  has 
been  employed  by  many  investigators  such  as  Abraham  and  Dare  (1971a,  b) , Burton 
(1970,  1972,  1973) > Dare  et  al  (1972),  Hale  and  Plummer  (197%),  Hoare  and  Pal 
(1971,  1975),  Kistenmacher  (197I+X  Lee  et  al  ( 1973a) > Nishioka  (1971),  and  Plummer 
and  Hale  ( 1972) . By  assuming  both  an  interaction  potential  between  atoms  and  a 
structural  model  for  the  cluster,  the  han  onic  oscillator  partition  function  of 
the  cluster  is  used  to  predict  the  thermodynamic  properties.  Although  this 
procedure  is  very  useful  for  studying  smAi 1 clusters  in  the  solid  phase,  it  is 
incapable  of  treating  liquids  because  of  the  simultaneous  assumption  of  inter- 
action potential  and  structure  configuration.  The  configuration  of  the  cluster 
cannot  be  given  independently  of  the  interaction  potential. 

• 

To  assure  this  consistency,  two  mathematical  methods  are  used:  the 
molecular-dynamics  method  (McGinty,  1973)  and  the  Monte  Carlo  method  (Lee  et 
al,  1973b) . The  molecular-dynamics  method  can  provide  both  static  and  kinetic 
properties  of  cluster  by  calculating  their  time- averaged  values.  The  Monte 
Carlo  method  gives  their  ensemble  average,  being  restricted  usually  to  static 
structural  properties.  The  partition  function  is  given  by  an  integral  of  con- 
figuration probabilities  with  respect  to  all  atoms  or  molecules  of  the  cluster 
system.  The  configuration  probability  is  a function  of  the  total  potential 
energy  of  the  interaction  between  the  atoms  of  molecules,  which  is  determined 
by  the  interatomic  or  intermole cular  potential.  Thus,  this  method  requires 
only  a knowledge  of  the  interatomic  or  intermolecular  potential. 

In  molecular  dynamics,  the  clusters  are  treated  as  molecules 
composed  of  spheres  with  assumed  intermolecular  potentials.  The  cluster 
partition  functions  are  then  evaluated  from  standard  statistical-mechanics 
relations,  using  approximations  such  as  the  harmonic  approximation  for  the 
vibrational  partition  functions  and  the  rigid  rotation  approximation  for 
the  rotational  partition  functions . Each  equilibrium  cluster  is  found  by 
relaxing  the  nonequilibrium  cluster  so  that  the  potential  energy  of  the  cluster 
takes  an  extremum.  From  the  partition  function,  thermodynamic  functions  for 
the  system  are  then  calculated,  and  the  equilibrium  concentration  of  clusters 
sure  obtained.  Such  molecular  dynamic  studies  of  microclusters  were  reported 
by  McGinty  (1973),  Briant  and  Burton  (1974,  1975a,  b,  1976a,  b),  Kristensen 
et  al  ( 1974),  and  Zurek  and  S chi  eve  (1978) . The  validity  of  such  a molecular 
approach  for  cluster  thermodynamics  was  studied  by  Courtney  et  al  (1974). 

In  the  Monte  Carlo  calculations,  a large  number  of  configurations 
are  generated  for  a molecular  cluster  system  in  which  the  clusters  are  defined 
by  a certain  criterion.  By  assuming  appropriate  intermolecular  potentials 
and  canonically  averaging  with  weighting  of  the  obtained  configurations 
corresponding  to  the  associated  potential  energy,  the  thermodynamic  properties 
of  the  cluster  are  calculated.  Such  calculations  were  performed  by  Lee  et  al 
(1973b),  Abraham  (1974a,  b) , Bonissent  and  Mutaftschiev  (1974) , Binder  (1975), 
Muller- Krumbhaan  and  Stoll  (1976),  Kaelberer  and  Etters  (1977),  and  Miyazaki 
et  al  (1977).  Detailed  considerations  of  the  definition  of  physical  clusters 
were  recently  given  by  Gill  (1977)  and  Gill  et  al  (1977) • 


Although  these  approaches  rest  on  the  conceptual  exactness  of  the 
statistical-mechanical  formulation  of  the  free  energy  of  cluster  formation, 
the  required  knowledge  of  the  intermolecular  potential  makes  it  difficult  to 
perform  the  requisite  calculations  for  arbitrary  substances.  Using  an  atomic 
model  of  clusters.  Burton,  and  B riant  (1977)  obtained  uni versed  correction 
factors , Iexact/- ci 3 which  were  independent  of  the  substance  studied.  The 
atoms  were  located  on  lattice  points  of  a face- centred  cubic  lattice  and 
interacted  only  through  nearest  neighbour  bonds.  Unlike  Abraham  ( 1969)  who 
divided  substances  into  the  two  categories  mentioned  above,  they  showed  that 
for  any  substances  at  low  supersaturation  and  low  temperature,  nucleation 
rates  should  be  lower  than  those  predicted  by  classical  theory,  while  at  high 
supersaturation  and  high  temperature,  higher  nucleation  rates  should  result. 

The  features  agree  with  observed  data.  However,  to  further  validate  the  theory, 
many  more  experimental  tests  will  be  required  over  a wide  range  of  supersatura- 
tion temperatures  and  liquids. 

Binary  and  multi- component  mixtures:  Multi- component  systems  consisting  of 
two  or  more  substances  of  condensable  vapour  are  very  important  for  chemical 
engineering  processes.  Nucleation  theories  of  the  one-component  system 
mentioned  above  can  be  extended  to  such  multi -component  systems.  Hirschf elder 
(1974)  extended  Reiss'  treatment  of  classical  theory  for  a binary  system  to  a 
multi -component  system.  Then  the  free  energy  of  formation  of  a cluster  of 
size  i is  given  by 

“i  ‘Xv*1!.  • “o  ) * "Si  <23> 

3 0 0 

where  |i0j  is  the  chemical  potential  per  molecule  of  the  j-th  component  in 
the  vapour  phase  having  the  same  composition  as  the  cluster,  pij  the  chemical 
potential  per  molecule  of  the  j-th  component  of  the  bulk  material  having  the 
same  phase  and  composition  as  the  cluster,  ij  the  number  of  molecules  of  the 
j-th  component  in  the  cluster  (i  = £ i.),  cr  the  surface  tension,  and  Si  the 

3 J 

surface  area  of  the  cluster.  In  binary  component  systems,  the  point  of 
unstable  equilibrium  for  the  critical  nucleus  corresponds  to  a saddle  point 
on  the  free-energy  surface  of  the  composition  plane. 

Recent  studies  in  this  field  were  made  by  Mirabel  and  Katz  (197*0, 
Heist  and  Reiss  (1974),  Shugard  et  al  (1974),  Hirschfelder  (1974),  Wilemski 
(1975a,  b),  Takahashi  et  al  (1975),  Reiss  et  al  (1976),  Bouland  (1977),  and 
Wegener  and  Wu  (1976) . These  theories  make  it  possible  to  determine  the 
value  of  the  nucleation  rate  in  relation  to  the  partial  pressures  of  the 
components.  For  binary  sulfuric-acid  water- vapour  systems,  even  for  humidi- 
ties less  than  100$,  nucleation  occurs  at  a low  sulfuric  acid  partial  pressure, 
because  the  mixture  can  be  supersaturated  with  respect  to  the  sulfuric  acid 
solution  (Heist  and  Reiss,  1974) . Using  raref action- wave  flows  of  H20  + C2H5QH 
binary  system  in  a shock  tube,  Wegener  and  Wu  (1976)  observed  a considerable 
reduction  in  the  critical  supersaturation  of  C2H5OH  with  increasing  concentra- 
tions of  H2O.  Wilemski  (1975b)  discussed  the  time  lag  needed  to  attain  steady- 
state  nucleation  in  binary  vapours.  Under  certain  conditions  where  the  concen- 
tration of  one  of  the  two  components  is  very  small,  the  time  lag  may  became  so 
large  that  the  classical  theory  of  nucleation  with  the  steady- state  assumption 
is  invalid.  Comparisons  between  these  theoretical  and  experimented,  studies  do 
not  always  show  good  agreement  due  to  the  uncertainty  of  the  vapour  pressure 
and  the  nature  of  the  resulting  chemical  compounds. 
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2.3  Heterogeneous  Nucleation 

Heterogeneous  nucleation  and  condensation  are  more  prevalent  phenomena 
than  homogeneous  nucleation  (Boucher,  1969,  Sigsbee,  1969)  in  the  cloud  physics 
as  well  as  in  metallurgical  fields.  However,  due  to  the  complexity  of  its 
physical  mechanisms  and  the  difficulty  of  mathematical  formulations,  the 
systematic  studies  have  made  little  progress  to  date.  The  main  sources  of 
heterogeneous  nuclei  in  the  atmosphere  are  smoke  and  vapour  from  fires  and 
various  industries,  dust  from  land  surfaces,  salt  from  oceans,  and  parti culate 
products  from  chemical  reactions.  The  chemical  and  physical  properties  of  these 
nuclei  are  extremely  varied  and  complicated. 

Heterogeneous  nucleation  and  condensation  can  be  characterized  as 
nucleation  under  the  existence  of  third  bodies,  B,  which  participate  in  the 
kinetic  process  as  follows: 


^ + B ^ ^ • B 

• B + ^ • B ?=*  Ag  • B + B 
• B *=*=  Ag  • B 
Ag  + B Ag  • B 

where  Ai  means  an  i-mer,  Aj/B  an  i-mer  strongly  combined  with  a third  body  B. 
Depending  on  the  phase  of  the  third  body,  heterogeneous  nucleation  and  conden- 
sation are  categorized  into  three  types:  B can  be  a (i)  solid,  (ii)  liquid, 
or  (iii)  gas  (molecule,  atom,  ion).  The  first  means  nucleation  from  the 
vapour  on  foreign  solid  particles.  The  monomer  or  i-mers  are  adsorbed  onto  or 
desorbed  from  the  solid  substrate.  In  the  second  case,  they  are  dissolved 
into  or  absorbed  into  the  liquid  material  depending  on  their  solubility. 

Hence  the  process  is  dominated  only  by  condensation  ( Saltancr  and  Tsirihin, 
1977) • The  third  case  usually  involves  nucleation  on  ions  and  has  been 
studied  in  some  detail. 

Nucleation  on  solid  particles:  The  process  of  heterogeneous  nucleation  on 
solid  particles  is  greatly  influenced  by  geometrical,  chemical  and  physical 
surface  properties  of  the  nuclei.  Taking  into  account  these  features,  several 
models  have  been  proposed  for  heterogeneous  nucleation  from  the  macroscopic 
point  of  view  (Hirth  and  Pound,  1963,  Sigsbee,  1969,  Boucher,  1969,  and 
Zimmels , 1976) . 


(24) 


Consider  monomers  of  a supersaturated  vapour  ircpinging  on  the 
surface  of  a solid  particle  of  radius  R,  to  make  a cap-shaped  cluster  embryo 
of  radius  r.  There  sure  two  ways  of  monomer  impingements  on  embryos:  a direct 
impingement  of  monomers  from  the  vapour  and  an  impingement  by  desorption  of 
adsorbed  monomers  (Fig.  3). 


As  for  the  homogeneous  case,  the  equilibrium  concentration  of 
embryos  of  size  i*  per  unit  surface  area  of  substrate  can  be  given  by 
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where  nsi  is  the  concentration  of  adsorbed  monomers  on  the  substrate  and 
&J*  is  the  free  energy  of  formation  of  an  embryo  of  size  i*  from  the  vapour 
phase  monomer  of  concentration  ns]_ . The  latter  concentration  can  usually  be 
expressed  in  terms  of  the  impinging  vapour  monomer  flux  1111%  and  the  meen 
monomer  stay  time  on  the  substrate  t as 


n , = a t 
si  a 


(26) 


where  aa  is  the  adsorption  coefficient  and  t is  given  by 


— = v exp 
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fd 

kl 


(27) 


where  is  the  desorption  free  energy  of  monomers  and  v is  the  vibrational 
frequency  of  the  adsorbed  monomer. 

In  the  direct  impingement  process  of  vapour  monomers,  the  impinge- 
ment rate  of  monomers  on  the  embryo  is  expressed  as  in  Eq.  (4), 


Csi  " \ 


nkT 


(28) 


where  ay  is  the  surface  area  of  the  emibryo.  The  nucleation  rate  per  unit  sur- 
face area  of  substrate  in  the  form  of  Eq.  (12)  is  then 


1 pv 


= z a*Vv 
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In  the  impingement  process  of  desorbed  monomers,  the  rate  at  which 
desorbed  monomers  impinge  upon  the  embryos  on  the  substrate  is  given  by 


= 1.1 


si 


t exp 
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(30) 


where  li  is  the  emibryo  perimeter,  is  the  mean  surface  diffusion  distance  and 
£GS  the  free  energy  of  activation  for  surface  diffusion.  The  nucleation  rate 
is  then 
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Using  reasonable  assumptions  that  £GS  « and  « ay,  both 

nucleation  rate  equations  can  be  expressed  in  the  same  form  as 
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where  n = 1 for  the  direct  impingement  of  monomers  from  the  vapour,  and  n ■ 2 
for  the  impingement  of  desorbed  monomers . 


From  the  macroscopic  surface-thermodynamic  point  of  view,  the  critical 
free  energy  of  the  embryo  formation  can  be  expressed  as  the  sum  of  contributions 
from  the  formation  of  the  liquid  bulk,  the  liquid- vapour  interface,  and  the 
solid-liquid  interface: 


“l  ■ “b  Ti  + V \ * (V»  - V)as 


(33) 


where  is  the  free-energy  difference  per  unit  volume  of  liquid  phase  from 
vapour  to  liquid  phase  being  equal  to  (nv  - oij  the  surface  free  energy 

of  the  interface  between  phases  i and  j,  vg  the  volume  of  liquid,  ay  the  surface 
area  of  liquid- vapour  interface,  as  the  surface  area  of  liquid- substrate  inter- 
face, and  subscripts  v,  i and  s refer  to  vapour , liquid  and  substrate,  respec- 
tively. From  geometrical  consideration,  vg,  ay  and  as  are  written  as 
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where  a,  by  and  bs  are  functions  of  r/R  and  the  contact  angle  0. 
By  using  the  usual  definition  of  the  contact  angle 
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and  substituting  Eqs.(34)  and  (35)  into  Eq.(33),  the  free  energy  of  formation  of  an 
embryo  of  radius  r,  on  a nucleus  of  radius  R,  can  be  expressed  as  a function  of 
r,  R and  cos0.  Differentiating  £Gi  with  respect  to  r to  evaluate  the  critical 
free  energy  gives  the  critical  radius  r*, 


r*  = - 


2j 


(36) 


The  free  energy  of  formation  of  a critical  embryo  is  then 
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where 
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Due  to  the  nature  of  the  factor  f(r/R,  cose) , it  is  seen  that  lower  values  of 
contact  angle  reduce  greatly  the  activation  energy  of  critical  embryo.  Finally, 
the  heterogeneous  nucleation  rate  per  unit  area  of  nucleus  substrate  can  be 
expressed  as 

Is'IoRH^yi;^{'>d(r- 1 )}“»{- ftt'"  - ^ 

where  IQ  means  the  homogeneous  nucleation  rate  per  unit  volume  of  vapour, 
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Nucleation  on  ions:  Clustering  and  condensation  of  vapour  atoms  or  molecules 
on  ions  have  been  investigated  since  Wilson  first  demonstrated  that  ions 
enhance  nucleation  in  cloud  condensation.  Ions  provide  a central  force  field, 
which  is  a simple  form  of  attractive  potential  resulting  in  the  formation  of 
heterogeneous  embryos.  Therefore,  a fundamental  understanding  of  ion- induced 
nucleation  would  provide  an  important  basis  for  further  developments  in  the 
theory  of  heterogeneous  nucleation.  In  an  attempt  to  explain  the  process  of 
nucleation  about  ions,  several  equations  have  been  proposed  by  White  and 
Kassner  (1971),  Mason  (1971),  Loeb  (1971),  Stachyorska  (1965),  Russel  (1969) > 
and  Castleman  and  Tang  (1972). 

Gaseous  Ions  attract  vapour  molecules  to  form  small  droplets  in  . 
equilibrium  with  the  vapour.  A larger,  single  charged  droplet  may  be  formed 
by  clustering  the  appropriate  number  of  monomers  to  one  droplet  of  radius  ra, 
containing  ia  molecules  (Fig.  4).  The  ion- field  energy  is  changed  due  to 
surrounding  dielectric  molecules . By  introducing  an  electrostatic  term  to 
account  for  the  charge,  the  macroscopic  free  energy  of  embryo  formation  is 
written  as 
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where  e is  the  liquid  dielectric  constant  and  e is  the  electronic  charge  in 
Gaussian  units.  At  r = ra,  = 0 and  <J£Gi/^r  = 0,  which  means  a met  as  table 
equilibrium.  The  first  two  terms  consist  of  a term  accounting  for  the  free 
energy  of  vaporization  of  the  condensing  molecules  and  a surface  term  re- 
presenting the  work  done  in  changing  the  surface  area  as  a result  of  the 
condensation  of  additional  molecules.  The  one  additional  term  accounts  for 
the  change  in  the  field  energy  about  an  ion  owing  to  the  condensation  of  a 
film  of  a given  radius  having  a dielectric  constant  e.  The  dielectric 
constant  of  the  surrounding  vapour  phase  is  generally  taken  to  be  close  to 
unity. 


The  critical  nucleus  ( ^Gi/br  = 0)  is  then  given  by 
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which  has  two  roots,  the  smaller  one  corresponding  to  ra  and  the  larger  to 
r*.  The  nucleation  rate  is  then 
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and  the  Zeldovich  factor  is  given  by 
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(43) 
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In  this  expression,  only  the  term  representing  different  ionic  radii 
accounts  for  different  ions,  and  it  does  not  always  predict  observed  variations 
in  the  nucleation  rate  when  the  sign  of  ions  is  changed.  White  and  Kassner 
(1971)  modified  the  term  to  account  for  dipole-dipole  and  dipole- quadrupole 
repulsions  between  the  clustering  molecules.  In  this  manner,  they  arrived  at 
a set  of  expressions  for  calculating  numerical  values  of  the  free-energy 
barrier  for  nucleation  about  both  negative  and  positive  ions.  The  values 
appeared  to  be  consistent  with  macroscopic  observations  of  the  relative 
saturation  ratios  required  for  nucleation  about  the  respective  ions  of 
different  sign . 

Russell  (1969)  extended  the  homogeneous  nucleation  theory  of  Lothe 
and  Pound  to  condensation  on  ions,  taking  the  clusters  as  free  particles  and 
accounting  for  the  corresponding  partition  functions.  The  free  energy  of 
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embryo  formation  is  then  expressed  as 


^ )/(  1 } <»> 

where  qtrans  and  q^t  are  molecular  translational  and  rotational  partition 
functions  of  the  entoryo  taken  as  a rigid  body,  and  qre_  is  a constant  term 
accounting  for  six  internal  degrees  of  freedom  replaced  by  the  translational 
.and  rotational  motions.  The  last  term  is  estimated  as  12kT  in(r/ra) , which 
leads  to  the  nucleation  rate 


(46) 


The  factor  (r */ra)1^  is  of  the  order  10^  in  most  cases.  Both  theories,  Eqs . 
(42)  and  (46),  are  compared  with  cloud  chamber  measurements  of  rain  limits, 
that  is,  the  super saturation  needed  to  produce  a rain  of  droplets.  Neither 
theory  predicted  the  different  catalytic  potencies  of  positive  and  negative 
ions. 


All  experimental  studies  of  ion-induced  nucleation  have  been  made 
in  cloud  chambers,  following  exposure  to  high-energy  X-rays  or  other  suitable 
ionization  sources,  and  measuring  only  macroscopic  properties  of  droplet 
formation.  Consequently,  the  molecular  aspects  of  the  phenomenon  and  espe- 
cially the  role  of  small  clusters  have  remained  obscure  and  many  conflicting 
statements  have  appeared  in  the  literature  concerning  the  actual  mechanisms 
involved*.  Mass-spectrometer  techniques  used  with  high-pressure  ion-sources 
have  yielded  a wealth  of  information  about  clustering  on  ions.  Water  has 
received  special  attention.  Clustering  rates,  equilibrium  distributions, 
and  other  thermodynamic  quantities  have  been  determined  for  small-size 
clusters  of  water  molecules  on  ions.  In  these  studies,  the  classical  liquid- 
drop  model  seems  woefully  unable  to  account  for  the  properties  of  these  i 
clusters,  although  it  may  be  more  realistic  as  clusters  grow  larger.  For 
example,  water  does  not  condense  around  such  ions  as  o£,  0+,  n£,  N*,  Ar+  or 
He+  as  had  been  proposed.  Instead,  these  ions  rapidly  react  with  water  to 
form  an  equilibrium  distribution  of  H+(H20)n,  which  serves  as  a nucleus  for 
further  condensation.  On  the  other  hand,  some  ions  such  as  Fb+,  Bi+,  Na+, 
and  Cl"  seem  to  form  centres  for  clustering.  Ion-molecule  reactions  are 
taking  place  during  the  course  of  the  usual  ion-induced  nucleation  experiment, 
generally  leading  to  the  formation  of  a complex  group  of  ions  of  unknown  com- 
position . 


To  account  for  these  features,  the  microcrystal  statistical  model 
developed  for  homogeneous  clusters  also  can  be  applied  to  vapour-ion  inter- 
actions. Casteman  and  Tang  (1972)  treated  the  clusters  as  harmonically-bound 
crystallities  and  calculated  their  parti tion  functions  directly  on  a high-speed 
cooputer . Once  the  partition  functions  are  determined  one  can  obtain  detailed 
thermodynamic  data  for  various  cluster  sizes  over  a wide  range  of  pressure  and 
a reasonable  range  of  tenperature.  With  the  molecular  dynamic  method,  B riant 
and  Burton  (1974)  calculated  the  thermodynamic  properties  of  the  clusters 
consisting  of  one  to  four  water  molecules  plus  one  ion  and  of  four  water 
molecules  plus  NH3,  assuming  an  effective  potential  for  water- water  and 


water-ion  interaction.  O'Brien  and  Robinson  (197*0  also  considered  two  sets 
of  clusters  of  spherical  shape  and  "grown"  clusters  of  an  ion  plus  different 
numbers  (l  ~ 86)  of  molecules.  Each  cluster  consists  of  Lennard- Jones  (6-12) 
spheres  (atcms  or  molecules)  with  a single  charge  ion  in  its  centre.  Mruisk 
et  al  (197**)  applied  the  Monte  Carlo  method  to  a system  of  water  vapour  plus 
Li+,  Na+,  K*,  Cl"  and  F"  ions  for  a cluster  of  one  to  six  molecules  of  H2O, 
and  obtained  sufficiently  good  agreement  with  experimental  data. 

2.4  Nonstationary  Nucleation 

In  the  case  of  condensation  occurring  in  a vapour  close  to  thermo- 
dynamic equilibrium,  the  features  of  the  condensation  process  sure  not  sensitive 
to  the  detailed  cluster- format! on  model  enployed,  and  they  are  described  by 
steady- state  kinetics.  However,  a steady  state  can  only  be  reached  when  all 
transient  effects  of  the  process  have  disappeared  from  the  system.  When  this 
relaxation  process  takes  a much  longer  time  than  the  characteristic  lifetime 
of  the  vapour  system  (as  for  a rapid  expansion  of  supersonic  free-molecular 
jets  or  in  the  associated  vapour  system)  the  transient  processes  play  an  impor- 
tant role  in  the  clustering  behaviour , For  this  purpose,  nonsteady  kinetics 
of  clusterization  must  be  considered. 

Earlier  approaches  to  this  problem  were  reviewed  by  Courtney  (1962). 
Employing  the  classical  liquid-drop  theory  with  a lower  limit  of  applicability 
as  a 19-molecule  cluster  and  assuming  a collision  frequency  for  the  rate 
constant,  Courtney  (1962)  solved  Eq.  (2)  and  displayed  the  time  dependence  of 
concentrations  of  various  size  clusters  (20  < i < 119) • Similar  studies  were 
done  by  Andres  and  Boundart  (1965)  and  Abraham  (1970,  1971)  for  11  < i < 110. 
The  latter  showed  an  overshoot  of  the  current  Ij.,  from  the  steady- sTate  current 
I*,  for  cluster  sizes  less  than  the  critical  size.  The  maximum  overshoot  is 
greater  in  magnitude  and  at  an  earlier  time  for  small  clusters.  Kortzeborn 
and  Abraham  (1973)  investigated  the  heterogeneous- nucleation  problem  of  water 
monomers  condensing  onto  a prescribed  nunfcer  of  ions  in  a similar  manner.  In 
this  case,  the  smallest  cluster  can  be  assumed  to  be  in  metastable  equilibrium 
with  a concentration  equal  to  that  of  the  ions.  The  time  dependence  of  cluster 
concentration  and  current  showed  the  same  functional  dependence  as  found  in 
the  homogeneous  case.  Firsch  and  Carlier  (1971) , Harstad  (1975) > Draine  and 
Salpeter  (1977)  and  Peak  (1978)  treated  the  system  of  cluster  reactions  as 
continuous  in  size  to  solve  Eq.  (2)  in  integral  form  and  calculated  the  relaxa- 
tion time,  accounting  for  the  effects  of  the  concentration  of  the  smallest 
cluster.  Shugard  and  Reiss  (1976)  applied  a stochastic  approach  to  the  non- 
steady nucleation  process.  The  growth  of  droplets  through  the  addition  of 
simple  molecules  can  be  analogous  to  a random  walk  of  walkers  on  a lattice  of 
discrete  size,  where  the  nunber  of  molecules  in  a drop  corresponds  to  the 
coordinate  of  the  lattice.  The  developed  theory  was  then  applied  to  the 
H2SO4-H2O  vapour  system,  which  had  quite  a long  relaxation  time  in  atmospheric 
conditions . 

As  mentioned  above,  in  free-molecular  Jet- expansions , the  kinetics 
of  formation  of  increasingly  larger  molecular  clusters  is  fundamentally 
necessary  to  characterize  the  condensation  process.  In  this  decade,  extensive 
efforts  have  been  made  in  this  area,  especially  for  dimer  formation.  This 
work  will  be  described  in  Section  4.4  in  connection  with  free-molecular  Jet 
flows.  Recently,  the  kinetics  of  clustering  was  studied  by  Pouring  (1977)  and 
Kanyukhov  and  Faizulaev  (1977)  considering  that  decay  and  growth  of  clusters 
follow  unimolecular  and  bimolecular  reactions. 
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Pouring  (1977)  applied  the  cluster  kinetic  theory  of  Buckle  (1969) 
to  the  condensation  of  water  vapour  in  air.  In  the  cluster  theory,  each  cluster 
is  an  assembly  of  atans  hound  by  London  forces.  The  probabilities  of  unlmole- 
cular  and  blmolecular  reactions  which  control  the  decay  and  growth  of  clusters 
are  determined  from  the  phase  integral  for  energy  configurations.  The  internal 
energy  in  a cluster  circulates  freely  between  all  degrees  of  freedom  at  some 
frequency.  When  the  energy  of  a surface-site  atom  exceeds  the  coupling  energy 
to  its  nearest  neighbours,  the  atom  is  lost  from  the  cluster.  On  the  other 
hand,  a cluster  can  grow  only  if  the  collision  energy  fails  to  dissociate  the 
incipiently-bound  atom  from  the  site  atom.  To  enploy  this  cluster  kinetics,  he 
assumed  molecular  parameters,  such  as  pair- interact! on  energy,  cluster  internal 
energy  circulation  frequency,  surface- site  to  internal-site  parameter , nearest- 
neighbour  coordination  numbers,  monomer  concentration,  and  cluster  size  of 
effective  numerical  Infinity.  After  an  initial  equilibrium  condition  is 
established,  the  resulting  relaxation  process  corresponding  to  changes  in  the 
thermodynamic  state  can  be  calculated  by  this  kinetic  theory.  It  was  found 
that  the  nonequilibrium  cluster  behaviour  is  dominated  by  the  behaviour  of  the 
internal  energy  redistribution  frequency. 

Zurek  and  S chi  eve  (1978)  used  molecular  dynamics  to  gain  insight 
into  the  details  of  cluster  formation.  Molecules  of  a two-dimensional  gas  of 
hard-sphere  well  discs  are  placed  initially  on  the  regular  lattice  with  equal 
speeds  and  random  directions  of  motion.  The  initial  condition  for  relaxation 
is  produced  by  a few  hundred  collisions.  Starting  from  this  initial  condition 
the  evolution  of  the  system  is  followed  by  Investigating  the  temperature 
dependence  of  the  mole  fractions  of  smai  i clusters  formed  in  an  adiabatically- 
isolated  system.  The  topological  aspect  of  clusters  also  was  studied.  The 
cluster  tends  to  have,  on  the  average,  a much  smaller  nuntoer  of  bonds  than  the 
most  compact  configuration. 


3-  CONDENSATION 

From  the  standpoint  of  nonsteady  kinetic  theory,  condensation  is  a 
successive  and  continuous  process  of  nucleation.  However,  in  usual  cases  of 
condensation,  the  nuclei  are  considered  to  consist  of  a sufficiently  large 
nunber  of  molecules  that  they  form  a liquid  droplet.  For  such  liquid  droplets, 
it  is  convenient  to  analyse  the  process  using  macroscopic  molecular  dynamics 
or  continuum  dynamics. 

3-1  Droplet  Growth 

In  the  process  of  droplet  growth,  the  molecular  motion  and  the  heat 
release  by  phase- transitions  must  be  taken  into  account  with  appropriate 
boundary  conditions  on  the  droplet  surface  for  incident  and  reflected  vapour 
and  gas  molecules. 

The  transfer  of  mass,  momentum  and  energy  to  droplets  depends  on 
the  Knuds en  maker  Kn  = Jfa/r^,  where  fa  is  the  mean-free-path  of  the  molecules 
in  the  gas  and  r±  the  characteristic  dimension  of  the  droplet.  The  mean-free- 
path  is  given  by 
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where  q is  the  dynamic  viscosity,  R the  gas  constant  and  p the  molecular 
weight.  The  limit  of  Kn  -»0  corresponds  to  the  continuum  regime  where  the 
transport  processes  are  governed  by  the  usual  fluid-dynamic  equations.  The 
mass  flux  of  the  vapour  to  a spherical  droplet  of  radius  r can  be  expressed 
as 

m = 4yr 2 — (C  - C ) 
o " r ' v«o  vry 


where  py  is  the  density  of  the  vapour,  D the  diffusion  coefficient  and  Cy,*, 
and  Cyr  the  concentrations  of  the  vapour  at  infinity  and  at  the  droplet 
surface,  respectively.  The  limit  of  Kn  -»  oo  corresponds  to  the  free-molecular 
regime  where  the  transfer  processes  are  controlled  by  the  interaction  of 
individual  gas  molecules  with  the  droplet.  The  mass  flux  of  the  vapour  to 
the  droplet  is  then  given  by  the  Hertz-Knudsen  relation 
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where  m is  the  molecular  mass,  py  the  vapour  pressure  at  the  droplet  surface, 
and  pa  the  saturation  vapour  pressure  at  the  droplet  temperature,  T^.  The 
first  term  in  the  parentheses  is  the  impingement  rate  of  vapour  molecule  on 
the  droplet  surface,  and  the  second  is  the  evaporation  rate. 


In  the  intermediate  regime,  these  expressions  must  be  modified  by 
an  amount  depending  on  the  deviation  from  the  foregoing  limits.  For  a case 
slightly  different  from  the  free-molecular  regime,  a sticking  coefficient,  ac, 
is  introduced,  and  the  droplet  growth  rate  can  be  expressed  as 


dr  m pv  ( pd  [t~  \ 

**  ’ c pi  jz* m V " PvV  Td  ) 
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where  pj  is  the  density  of  liquid  of  the  droplet.  By  introducing  an  effective 
condensation  coefficient 


(48) 


Eq.  (47)  can  be  reduced  to 


dt  c pi  4RSS 


(47a) 


The  sticking  coefficient  ctc  is  the  ratio  of  the  number  of  vapour  molecules  (or 
atoms)  that  become  part  of  the  liquid  phase  to  the  toted,  number  of  molecules 
that  impinge  on  the  surface  or  the  fraction  of  incident  vapour-phase  molecules 
which  are  trapped  in  the  liquid- vapour  surface  potential-well  after  a single 
collision  with  the  surface.  The  condensation  coefficient  is  synonymous  with 
the  sticking  coefficient.  This  concept  and  the  concept  of  thermal  accommodation 
coefficient  form  the  basis  of  the  equations  for  droplet  growth  in  this  regime. 
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Analytic  expressions  for  droplet  growth  consider  the  diffusion  of  vapour  toward 
the  droplet  and  the  conduction  of  heat  away  from  the  droplet.  The  sticking 
coefficient  ac  in  the  equation  modifies  the  binary  diffusion  coefficient  for 
the  vapour  diffusing  through  the  noncondensing  gas.  The  thermal  accommodation 
coefficient  at  modifies  the  thermal  conductivity  of  the  noncondensing  gas.  The 
latter  coefficient  is  a measure  of  the  efficiency  with  which  the  noncondensing 
gas  removes  heat  from  the  droplet  surface.  In  water-droplet  growth  as  applied 
to  the  atmosphere,  it  appears  reasonable  to  use  an  accommodation  coefficient 
of  unity.  However,  the  sticking  coefficient  varies  and  depends  on  a number  of 
factors  such  as  the  degree  of  contamination  of  the  vapour  and  the  nature  of  the 
noncondensing  gas. 

Values  of  the  sticking  coefficient  can  be  obtained  theoreti cally 
from  kinetic  theory  by  solving  the  Boltzmann  equation  in  the  vapour  phase 
subject  to  the  appropriate  boundary  conditions  at  the  surface  of  the  droplet. 

The  solution  can  be  expanded  in  terms  of  the  inverse  Knudsen  number.  The 
coefficients  of  this  expansion  are  determined  by  collision  integrals.  This 
problem  was  studied  by  Kelly  and  Sengens  (1972,  197*+),  and  Loylka  (1973)  with 
a linearized  BGK  model. 

Values  of  the  sticking  coefficient  have  also  been  obtained  experi- 
mentally by  several  investigators  for  a variety  of  condensation  and  evaporation 
phenomena  and  experimental  techniques.  According  to  Mills  and  Seban  (1967), 
the  value  ranged  from  0.006  to  unity.  By  means  of  thermal  diffusion  chambers, 
Chodes  et  al  (197*0,  Gollub  et  al  (197*+)  and  Sinnarwall  et  al  (1975)  obtained 
an  average  value  of  0.036.  Vietti  and  Fastook  (1976)  measured  water-droplet 
growth  in  air  and  CO2  atmospheres  using  an  expansion  cloud  chamber  with  a light 
scattering  technique  for  saturation  ratios  s = 1.5  ~ 3*5*  Comparing  with  the 
Fukuta-Walter  theory  (Eq.  53) > they  obtained  a value  of  the  sticking  coefficient 
ac  *»  1 for  air  and  ac  = 0.1  ~ 1 for  a CO2  atmosphere  (r  < 2 pm) . For  larger 
droplets  in  C02,  values  smaller  than  0.036  were  observed.  With  a diffusion 
cloud  chamber  and  a light  scattering  method,  Gollub  et  al  (197*+)  obtained 
values  for  the  condensation  coefficient  and  the  thermal  accommodation  coeffi- 
cient for  small  saturation  ratios  (s  = 1.02  ~ 1.05)  as  at  = 1 and  ac  = 0.01  ~ 0.1 
for  the  3 ~ 7 pm  radius  range. 


From  continuum  fluid- dynamics,  several  expressions  are  proposed 
for  the  droplet  growth  in  the  intermediate  regime.  When  the  droplet  population 
is  sufficiently  low,  the  droplet  and  vapour  system  can  be  treated  as  a single 
droplet  surrounded  by  an  infinite  vapour  and  noncondensible  gas  mixture. 
Gyarmathy  (1963)  considered  mass  and  heat  transfer  about  a liquid  drcplet  of 
radius  r in  a mixture  of  t vapour  and  noncondensible  gas,  which  lead  to 
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where  m is  the  condensation  rate  of  vapour  mass  per  unit  time,  being  equal  to 
d(*+Tr/3r3p^)/dt,  Tr  the  surface  tenperature-  of  the  droplet,  T the  bulk  gas 
temperature,  pyr  and  py  the  vapour  pressure  at  the  droplet  surface  and  at  ' . 
infinity,  respectively,  psr  the  saturation  vapour  pressure  at  Tr,  py.  the 
density  of  the  vapour,  Rv  the  gas  constant  of  the  vapour,  cr  the  surface  tension, 
L the  latent  heat  of  vaporization,  and  Jt  and  (c  the  transfer  coefficients  of 
heat  and  mass  modified  for  the  Knudsen  flow: 
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where  X is  the  heat  conductivity,  D the  binary  diffusion  coefficient,  p the 
ambient  pressure,  pg  the  pressure  of  the  gas,  and  a^  and  ac  are  given  by 
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where  a^  and  a^v  are  the  thermal  accommodation  coefficients  for  the  mixture 
and  the  vapour , respectively,  tj  the  dynamic  viscosity,  Cp  the  specific  heat, 
7 the  specific  heat  ratio,  p the  density  of  the  mixture  and  p and  uy  the 
molecular  weights  of  the  mixture  and  the  vapour,  respectively.  From  these 
relations  and  the  Clasius-Clapeyron  relation 
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the  droplet  growth  can  be  expressed  as 
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where  rff  = 2<r/(^vRvT)  (in(py/pa)  }_1.  Enploying  this  droplet  growth  law  with 
the  classical  theory  of  homogeneous  nudeation  rate,  Eq.  (l8),  Barschdorff 
et  al  (1976)  showed  good  agreement  between  the  theoretical  and  experimental 
results  in  nozzle  and  shock-tube  flows  with  condensation. 

Fukuta  and  Walter  (1970)  derived  a similar  expression  for  the 
droplet  growth  law  in  the  case  of  py,,,  = 0 (pg«,  = p) , in  which  atv  in  *c  is 
replaced  by  the  condensation  coefficient  ac 
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Measuring  the  growth  rate  of  water  droplets  in  air,  argon  and 
helium  atmospheres,  in  an  expansion  cloud  chamber  with  a light  scattering 
method,  and  comparing  the  results  with  various  expressions  of  droplet  growth, 
Vietti  and  Schuster  (1973a,  b)  were  not  able  to  find  a best  fit  for  their 
data. 

For  densely  populated  droplets,  the  interactions  between  droplets 
must  be  taken  into  account  due  to  the  Brownian  motion,  electric  charge  and 
turbulent  flow.  In  highly  dispersed  systems,  coagulation  due  to  the  Brownian 
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motion  as  well  as  the  environmental  variation  due  to  vapour  depletion  and 
heat  release  of  condensation  will  predominate.  On  the  other  hand,  turbulent 
flows  will  play  an  important  role  as  the  inter-droplet  distances  increase. 
The  last  effect  was  studied  as  an  extension  of  a single  droplet  growth  by 
Wagner  and  Pdhl  (1975)  and  Wagner  (1975) • 

3.2  Condensate  Mass 


To  calculate  the  flow  variables,  it  is  necessary  to  know  the  overall 
condensate  mass  as  a function  of  time  and  space.  For  the  nonsteady  process 
of  nucleation  with  a large  relaxation  time,  the  history  of  the  condensate  mass 
must  be  conside^jsd  vising  nonsteady  kinetic  theory.  For  most  systems  with 
fana.ll  relaxation  times,  it  can  be  considered  that  liquid  droplets  with  the 
critical  radius  are  spontaneously  nucleated  corresponding  to  the  vapour  state 
and  then  grow  to  larger  droplets  obeying  a growth  law.  Through  vapour  conden- 
sation, droplets  nucleated  at  a time  t with  critical  radius  r*(  t)  grow  to 
droplets  of  radius  r(t)  at  a time  t.  Since,  during  the  time  interval  t to 
t + &t,  the  number  of  nucleated  droplets  is  l(-r)dT,  the  overall  condensate 
mass  per  unit  volume  of  the  mixture  at  time  t is  written  as 

t 

M = J ^ Pj|  I(r)r(t,  T)3dT  (55) 
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The  rate  of  production  of  the  overall  condensate  mass  is  then 

t 

« = 31  PiKt)r*(t)3  + J pf(r)  ^ (r(t,  t)  }3dr  (56) 
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When  nucleation  and  condensation  proceed  heterogeneously  on  the 
surface  of  foreign  solid  or  liquid  particles,  an  additional  contribution  must 
be  added  to  the  above  relation.  In  order  to  describe  the  nunfcer  deisity  and 
size  distribution  of  such  foreign  nuclei,  a distribution  function  with  respect 
to  the  characteristic  size  of  particles  of  an  equivalent  radius  R,  can  be 
introduced.  The  distribution  function  N(R)  multiplied  by  dR  characterizes  the 
nusfcer  of  particles  lying  in  the  size  range  R to  R + dR  per  unit  volume  of 
vapour, 
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A nucleus,  particle  has  the  surface  area  S0  at  the  initiation  of  nucleation. 
At  time  t,  the  area  comes  to  be  wetted  by  condensed  droplets,  say,  Ss»and  the 
surface  area  of  a nucleus  particle  on  which  embryos  can  nucleate  is  then 
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When  the  nuntoer-distribution  function  of  foreign  nuclei  is  N(R)  and  the 
nucleation  rate  per  unit  surface  area  of  a nucleus  particle  is  Is,  the 
total  nucleation  rate  per  unit  volume  of  vapour  is  given  by 
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The  wetted  surface  area  Ss  is  given  by 
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where  ag(t,  t)  is  the  interfacial  area  between  the  nucleus  particle  and  the 
condensate  droplet  which  is  formed  at  time  t (<  t) . The  mass  produced  during 
the  time  r and  t + dT  is 
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and  grows  at  time  t up  to 
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where  vg  represents  the  volume  of  the  condensate.  Thus,  the  overall  conden- 
sate mass  per  unit  vapour  volume  is  given  by 
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and  the  rate  of  production  of  the  liquid  phase  by 
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When  condensation  of  the  vapour  takes  place  on  the  whole  surface  area  of  the 
foreign  particle  as  in  the  case  of  liquid  particles,  it  can  be  assumed  that 
Ss  = 0 and  SQ  ■ byr2. 

Combining  Eqs.  (58)  and  (62)  yields  the  final  expression  for  the 
rate  of  production  of  the  liquid  phase  as 

M = I(t)  ^ pjr*3  +J  I(t)  (3*  Pf3  ) dT 
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4 . FLOWS  WITH  NUCLEATION  AND  CONDENSATION 


For  a given  initial  condition  of  saturated  or  supersaturated  vapour, 
a finite  time  is  required  for  the  molecules  to  achieve  appreciable  clusterisa- 
tion or  critical  nucleation.  Such  relaxation  times  of  nucleation  are  of  the 
order  of  microseconds  (Andres,  1969)  • Thus,  the  vapour  flow  proceeds  without 
condensation  beyond  the  saturation  state,  being  supercooled  to  a certain  extent 
before  the  occurrence  of  observable  nucleation  (Fig.  5)  • In  cloud  chambers, 
supersonic  nozzles,  and  shock  tubes,  the  characteristic  time  of  flow  in  which 
an  appreciable  change  of  the  thermodynamic  state  occurs  is  of  the  order  of 
milliseconds.  The  cooling  rates  of  these  flows  usually  obtained  are  shown  in 
Table  1.  For  flows  where  the  thermodynamic  state  changes  sufficiently  slowly 
compared  with  the  kinetic  process  of  clusterization,  the  kinetically  steady 
state  can  be  established  after  a very  short  induction  time  for  nucleation, 
and  the  flow  is  described  by  a quasi- steady  theory  of  nucleation.  In  such 
cases,  nucleation  and  condensation  can  be  regarded  only  as  a heat- release 
process  and  the  physical  properties  of  the  flow  can  be  determined  from  well- 
established  methods  of  gasdynamics.  Consequently,  such  flows  have  been  used 
for  testing  and  understanding  nucleation  and  condensation  theories  and  to  provide 
macroscopic  theories  of  nucleation  sufficiently  useful  for  practical  problems  of 
condensation . 

Under  rapid  expansions  of  free-molecular  beams  produced  by  expansion 
through  an  orifice  or  nozzle  into  high  vacuum,  the  quasi- steady  assumption  of 
clusterization  is  no  longer  valid  and  a kinetic  treatment  must  be  used.  In 
turn,  the  kinetic  processes  can  be  studied  in  such  a flow  system,  and  detailed 
information  can  be  obtained  on  small  clusters. 

4.1  Cloud  Chambers 

— 

(i)  Expansion  Cloud-Chamber: 


The  cloud  chamber  commonly  discussed  by  many  investigators  is  the 
Wilson  expansion  chamber.  Here,  a mixture  of  vapour  and  inert  carrier  gas  is 
expanded  rapidly  from  a known  initial  volume  to  a known  final  volume  under 
well- controlled  conditions.  The  rapid  expansion  is  usually  obtained  by 
actuating  a piston  or  by  opening  a quick-acting  valve  to  a lower  pressure. 

The  final  state  in  the  chamber  is  chosen  to  be  supersaturated.  By  increasing 
the  final  supersaturation,  the  critical  limit  of  supersaturation  is  determined 
beyond  which  the  supersaturated  state  can  no  longer  be  sustained  owing  to 
onset  of  condensation.  At  this  point,  the  nucleation  rate  per  cubic  centi- 
meter per  second  is  considered  to  lie  in  the  neighbourhood  of  I = 1 drop/(cm3s) 


For  the  study  of  homogeneous  nucleation,  the  chamber  is  cleared  of 
heterogeneous  condensation  nuclei  by  repeated  expansion  with  the  same  mixture 
and  removing  the  drops  condensed  on  these  nuclei.  Possible  condensation  of 
vapour  on  ions  is  avoided  by  enploying  a suitable  electric  field  and  shielding 
the  chamber  against  low-energy  radiation.  The  critical  supersaturation 
required  for  homogeneous  nucleation  is  calculated  by  assuming  that  nucleation 
occurs  at  the  maximum  supersaturation  achieved  through  expansion  and  assigning 
An  appropriate  duration  time  to  this  state  - denoted  as  the  sensitive  time. 
Super saturation,  however,  is  continuously  and  rapidly  varying  with  time,  and 
additional  nucleation  and  droplet  growth  are  simultaneously  and  continuously 
occurring.  The  analysis  of  the  results  obtained  by  photographing  the  droplets 


or  from  light  scattering  techniques  must  be  studied  using  a nucleation-rate 
law  and  a droplet-growth  law.  A cloud  chamber  of  thiB  type  for  the  study  of 
homogeneous  nucleation  was  developed  to  a high  degree  of  refinement  by  Allard 
and  Kassner  (1965),  and  Kassner  and  Schmitt  (1966).  They  used  a programmed 
expansion  - a large  fast  expansion  followed  by  a small  slow  expansion  and 
a small  fast  compression.  The  programmed  expansion  produces  a roughly  para- 
bolic pulse  of  supersaturation.  The  small  slow  expansion  offsets  the  increase 
in  the  temperature  (decrease  in  supersaturation)  of  the  interior  gas  of  the 
chanter  which  is  brought  on  by  the  heated  boundary  layers  of  the  gas  near  the 
walls  and/or  around  growing  drqplets.  The  small  compression  reduces  the 
supersaturation  to  an  extent  so  as  to  suppress  the  additional  nucleation. 

At  the  same  time,  the  super saturation  is  kept  high  enough  to  insure  rapid 
droplet  growth.  The  pulse  technique  can  provide  sampling  times  which  are 
considerably  shorter  than  the  natural  "sensitive  time"  of  the  chamber.  With 
photographic  records  of  the  total  droplet  population, which  are  taken  after 
the  droplets  have  grown  to  visible  size,  the  relationship  between  nucleation 
rate  and  supersaturation  is  estimated  in  the  following  way.  First,  an 
espirical  nucleation-rate  versus  supersaturation  curve,  l(s),  is  assumed. 

The  curve  is  integrated  over  the  pulse  shape  of  the  experimented,  supersaturation 
to  predict  the  droplet  population 


u 

n(t)  = J I(s)dt 


The  empirical  curve  l(s)  is  adjusted  until  the  predicted  value  of  the 
population  agrees  with  the  one  obtained  photographiceLlly.  The  fined, 
nucleation-rate  versus  supersaturation  curve  includes  the  effect  of  vapour 
depletion  due  to  nucleation  and  droplet  growth,  which  may  not  necessarily  be 
considered  in  the  nucleation  theory.  To  avoid  the  latter  effect,  the  sampling 
should  be  taken  in  an  infinitesimally  short  time  or  the  droplet  growth  by 
diffusion  should  be  taken  into  account  under  the  time-varying  super saturation 
owing  to  vapour  depletion  as  well  as  continuous  expansion  or  compression . 

Such  an  attenpt  was  made  by  Lund  and  Rivers  (1966) . They  modified 
a droplet- growth  law  in  cloud  physics  used  by  Fletcher  and  others  for  the  growth 
and  evaporation  of  cloud- chamber  droplets.  Using  the  classical  nucleation 
theory  and  the  droplet  growth  law,  following  the  step-by-step  growth  of 
drqplets  and  correcting  nominal  supersaturations  obtained  from  pressure  records 
for  vapour  depletion,  they  predicted  droplet-population  and  droplet- size  in 
reasonable  agreement  with  observations.  However,  an  exact  treatment  of  this 
problem  is  complicated  by  the  necessity  of  solving  the  nonsteady  equations  of 
heat  and  mass  transfer  with  moving  boundary  conditions  which  are  time -dependent. 
For  this  purpose,  the  technique  vised  for  a shock- tube  raref action- wave  flow  can 
be  enployed  (see  Section  4.3) 

In  addition  to  such  nonstationary  processes,  the  results  obtained 
from  an  expansion  cloud- chamber  depends  on  the  degree  of  adiabatic! ty  of  the 
apparatus,  size  of  the  chamber,  "sensitive  time"  of  esqpansion,  and  observation 
techniques.  These  ccnplexities  have  been  circumvented  by  the  development  of 
the  diffusion  cloud- chamber . 

(ii)  Diffusion  Cloud- Chamber 

During  the  past  decade,  Katz  and  his  coworkers  (1967,  1970a,  1975, 
1976)  have  refined  the  thermal-diffusion  cloud- chamber  to  such  an  extent  that 


the  chamber  is  now  one  of  the  best  quantitative  methods  for  studying  the 
nucleation  process  in  a gas.  As  shown  in  Fig.  6,  it  consists  of  a warm  bottom 
plate  and  a cooler  top  plate.  A pool'  of  working  fluid  under  investigation 
covers  the  lower  plate.  The  space  above  the  pool  is  filled  with  a light, 
inert  gas  such  as  helium.  The  side  walls  of  the  chamber  form  a cylinder  and 
can  be  constructed  of  glass,  quartz  or  other  materials  containing  windows 
through  which  optical  observations  can  be  made.  Vapour  evaporates  free  the 
surface  of  the  pool  and  diffuses  to  the  upper  plate  where  it  condenses  at 
lower  temperatures.  The  upper  plate  is  slightly  beveled  and  the  condensate 
flows  to  the  side  walls  where  it  returns  to  the  lower  pool.  The  total 
density  of  the  vapour- carrier-gas  mixture  decreases  from  the  bottom  to  the  top 
to  prevent  natural  convection.  Since  the  saturation  vapour  pressure  of  the 
mixture  varies  exponentially  with  tenperature  and  the  partial  pressure  of  the 
vapour  is  nearly  linear,  a supersaturated  state  can  be  obtained  through  the 
chamber.  At  given  tenperatures  of  the  upper  and  lower  plates,  a steady  state 
is  set  up  in  which  supersaturation  passes  through  a maximum  at  a point 
relatively  distant  from  either  plate.  In  order  that  the  distribution  of 
tenperature  and  supersaturation  in  the  chamber  should  be  easily  predicted, 
the  chamber  must  be  designed  to  produce,  to  a high  degree  of  approximation, 
one- dimensional  diffusion  of  heat  and  mass  and  it  must  always  operate  under 
conditions  in  which  very  few  droplets  form. 

A single  experiment  consists  of  setting  the  tenperature  gradient 
such  that  the  chamber  is  visually  observed  to  be  nucleating  at  a rate  of  about 
2 ~ 3 drops  cm- 3 sec“l  and  then  measuring  the  total  pressure  and  the  tenpera- 
tures at  the  surface  of  the  evaporating  pool  and  at  the  top  plate.  The 
tenperatures  of  both  plates  are  then  decreased  or  increased  and  the  experi- 
ment is  repeated  in  such  a way  that  the  same  rate  of  nucleation  is  always 
observed.  The  critical  supersaturation  at  a given  tenperature,  the  super- 
saturation which  just  initiates  nucleation  at  the  assigned  rate  (of  1 ~ 10  droplets 
cm"3  sec-1)  can  be  obtained  as  the  envelope  of  the  individual  supersaturation 
curve  plotted  versus  the  tenperature  for  each  experiment  (Fig.  7). 

Using  this  technique,  Katz  and  his  coworkers  carried  out  an 
extensive  study  to  determine  the  critical  supersaturations  for  homogeneous 
nucleation  of  various  substances  such  as  ethanol,  methanol,  n- hexane  and 
water  in  the  tenperature  range  235  to  295  K (1967),  n-hexane,  n-heptane,  n- 
octane  and  n-nonane  in  the  tenperature  range  225  to  330  K (1970a),  benzene, 
toluene,  o- xylene  and  n-butylbenzene  over  as  much  as  a 150  K- range  of  tenper- 
atures (1975) > and  carbon  tetrachloride  (CCI4),  chloroform  (CHCI3) , Freon  11 
(CCI3F)  and  1-1-2-2-tetrachloro-ethane  (C2H2C1J|)  over  as  much  as  a 90  K- range 
of  tenperatures  (1976) . The  results  obtained  were  compared  with  predictions 
from  homogeneous  nucleation  theories.  The  classical  theory,  Eq.  (l8) , was 
found  to  be  in  excellent  agreement  without  any  correction  for  the  dependence 
of  surface  tension  on  radius  even  at  radii  as  small  as  9-7  A containing  l6 
molecules  of  n-butylbenzene  at  215  K or  11.7  A containing  2k  molecules  of 
n-nonane  at  230  K.  However,  in  shock  tube  flows  where  the  particle  radius 
is  less  than  10  A,  it  was  necessary  to  "adjust"  the  values  of  surface  tension 
for  a good  fit  of  analysis  with  experiment  (Sislian,  1976  and  Sislian  and 
Glass,  1977) • The  results  flw  chloroform  and  Freon  11  obtained  with  a 
supersonic  nozzle  showed  a good  agreement  with  the  Lothe-Pound  theory,  Eq.  (22) 
(Jaeger  et  al,  1969,  Dawson  et  al,  1969) . However,  now  the  results  from  a 
diffusion  chamber  were  found  to  be  in  very  good  agreement  with  classical 
theory,  Eq.  (l8) . By  using  hydrogen,  helium,  neon,  nitrogen  and  argon  as  the 
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carrier  gas,  the  effects  of  quantities  and  species  of  carrier  gas  were  also 
studied.  As  long  as  the  total  pressure  was  greater  than  about  2.5  times  the 
equilibrium  vapour  pressure  at  the  top  of  the  bottom  plate,  the  amount  of 
carrier  gas  was  found  to  have  no  effect  on  the  observed  critical  supersatura- 
tion. 

The  experimental  data  for  water  vapour  obtained  by  Katz  and  Ostermier 
(1967)  exhibited  considerable  scatter  by  comparison  with  their  other  results. 
Keist  and  Reiss  (1973)  considered  that  this  arose  because  the  water  which 
condensed  on  the  top  plate  did  not  wet  the  metal  surface  but  formed  large  drops 
which  fell  through  the  chamber  to  make  disturbances.  By  covering  the  top  plate 
with  a thin  layer  of  Pyrex  glass  wool  so  that  water  wetted  the  glass  wool  to 
drain  to  the  side  walls,  they  succeeded  in  eliminating  the  scatter  and  obtained 
reliable  data  for  homogeneous  nucleation  of  water.  In  this  case,  the  classical 
nucleati on-theory  provided  excellent  agreement  with  experimental  results . It 
should  be  noted  that,  by  contrast  with  nozzle  and  shock- tube  flows,  almost 
all  the  results  from  cloul  chamber  are  in  good  agreement  with  the  classical 
droplet-model  theory.  At  present,  it  is  not  clear  whether  this  is  due  to  the 
physical  features  of  cloud  chamber  or  to  the  limited  experimental  conditions. 


The  diffusion  cloud- chamber  has  the  following  advantages:  (a)  it 
is  a steady- state  device  such  that  a stationary  state  of  a specified  nucleation 
rate  can  be  used  to  determine  the  critical  super saturation;  (b)  the  chamber  is 
self- cleaning.  Once  e dust  particle  has  nucleated  a droplet,  it  rapidly  grows 
to  visible  size  and  Tills  into  the  liquid  pool;  (c)  the  chamber  effectively 
has  no  walls;  (d)  nucleation  on  ions  or  those  caused  by  the  passage  of  an 
energetic  particle  are  readily  distinguishable  from  homogeneous  nucleation 
However,  its  applicability  is  limited  by  the  following  factors:  (i)  the  sub- 
stance to  be  studied  must  be  liquid,  with  the  temperature  below  the  critical 
point  and  above  the  triple  point,  at  both  the  bottom  and  top  plates;  (ii) 
the  heat  and  mass  transfer  equations  require  an  accurate  knowledge  of  the 
thermal  conductivity  and  the  diffusion  coefficient  of  the  vapour- carrier-gas 
mixture  as  a function  of  both  temperature  and  composition. 

The  diffusion  cloud  chamber,  however,  because  of  its  advantages, 
can  be  used  for  other  studies  such  as  nucleation  on  ions,  in  binary  or  multi- 
component systems.  Furthermore,  the  nucleation  and  growth  processes  have 
the  possibilities  of  being  applied  to  the  detection  and  amplification  in  cheo- 
ical-physical  phenomena  involving  neutral  molecules.  Reiss  et  al  (1977)  have 
uncovered  many  experimental  difficulties  in  adapting  the  diffusion  chamber  for 
this  purpose  but  nevertheless  obtained  some  quantitative  results. 

4 .2  Nozzle  Flows 

Nucleation  and  condensation  in  nozzle  flows  are  important  in 
connection  with  the  design  and  performance  of  supersonic  and  hypersonic  wind 
tunnels,  rocket  nozzles,  steam  turbines,  and  cooling  towers.  Extensive  theore- 
tical and  experimental  studies  have  been  reported  since  Stodola  in  1927 . As 
a result,  ncezle  condensation  is  now  understood  in  same  detail.  Recently,  the 
nozzle-flow  technique  was  applied  successfully  to  studies  of  homogeneous - 
nucleation  and  condensation  problems.  Some  favourable  features  of  nozzle-flow 
experiments  are  that  (a)  the  state  is  steady  so  that  the  entire  history  of 
expansion  is  readily  obtained,  (b)  the  flow  between  the  wall  boundary  layers 
may  generally  be  considered  as  a uniform  adiabatic  flow,  (c)  nucleation  is 
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unlikely  to  occur  in  the  boundary  layers  since  these  are  largely  superheated 
due  to  friction  effects,  (d)  the  rate  of  condensation  can  be  estimated  by 
measuring  wall  static  pressures.  A good  review  of  the  subject  was  given 
recently  by  Wegener  (1969,  1975 , 1977),  where  relevant  previous  papers  up 
to  1976  are  listed. 

When  a vapour  or  vapour  mixture  flows  along  the  axis  of  a super- 
sonic nozzle,  the  vapour  pressure  decreases  with  decreasing  temperature  more 
rapidly  than  the  expansion  pressure  down  to  or  below  the  saturation  pressure. 
Further  expansion  may  lead  to  condensation,  which  causes  first  a slight 
deviation  of  static  pressure  from  the  isentropic  eaqpansion  and  then  an  almost 
abrvpt  increase  in  the  pres stare  in  a short  distance  (Fig.  8) . In  this  conden- 
sation zone,  the  maximum  heat  release  is  attained.  Further  expansion  leads 
to  two-phase  flow  of  an  asymptotically  saturated  vapour  and  liquid  droplets. 
Depending  on  prevailing  stagnation  conditions,  nucleation  and  condensation 
can  also  occur  in  subsonic  nozzles  (Barschdorff  et  al,  1976) . 

Typical  nozzles  used  for  condensation  studies  (except  free-molecular 
beams)  have  high  Reynolds  nunfcers  lCr  to  10°.  Therefore,  viscous  effects  are 
negligible  and  the  boundary  layers  are  confined  to  a very  thin  region  close 
to  the  walls.  In  'addition,  owing  to  negligible  wall  curvature  one  can  obtain 
a uniform  velocity  distribution  across  the  cross-sections.  For  such  flows  at 
high  Reynolds  numbers  and  low  supersonic  Mach  .nuntoers,  the  concepts  of  one- 
dimensional inviscid  continuum  flow  are  applicable.  Employing  this  assumption, 
Hill  (1966)  showed  that  the  features  of  condensation  flows  in  a nozzle  were 
well-predicted  by  the  gasdynamic  equations  when  coupled  with  the  classical 
nucleation  theory. 


For  a mixture  of  a condensible  vapour,  its  liquid  condensate,  a 
noncondensible  inert  gas  and  foreign  liquid  or  solid  particles  (suffixes  v, 
l , g and  p) , the  equations  of  mass,  momentum  and  energy  are: 

h + V “ 0 (64) 
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with  w - mass-flow  rate,  u ~ one- dimensional  axial  flow  velocity,  p - pressure, 
h - enthalpy  and  A - cross-sectional  area  of  the  nozzle.  The  cross-sectional 
area  is  adjusted  by  a small  enpirical  correction  for  the  boundary- layer  effect. 


When  the  no- slip  condition  of  velocity  and  temperature  between  the  gas, 
the  condensate  and  the  particle  applies,  the  equations  can  be  reduced  to  simple 
inviscid  one- dimensional  form: 
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where  the  density  p and  the  specific  heat  CpQ  are  defined  by 
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with  pi  and  pi  as  the  mass  of  the  condensate  (=M)  and  the  foreign  particles 
per  unit  volume  of  the  mixture,  respectively;  cp  - specific  heat  at  constant 
pressure,  L - latent  heat  of  vaporization  of  the  liquid,  and 
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w is  the  mass  fraction  of  the  condensable  vapour  and  the  condensate  - a 
constant  throughout  the  expansion;  g and  gp  are  the  mass  fraction  of  the 
condensate  and  the  particles,  respectively.  The  presence  of  the  term  gL  in 
the  energy  equation  means  that  vapour  condensation  is  equivalent  to  a heat 
addition  per  unit  mass  of  the  flow. 

The  equation  of  state  can  be  written  as 


where  p is  the  molar  weight.  For  certain  fluids,  say  high-pressure  steam 
or  metal  vapour , the  equation  of  state  may  deviate  from  an  ideal  gas  and 
must  be  expressed  in  a more  complex  form. 

The  mass  fraction  of  condensate  is  given  by  Eq.  (55) , 
x 

g(x)  =J  31  P/i(t)r(x,e)3  (?1) 

-00 

For  the  nucleation  rate  l( {) , Eq.  (l8)  can  be  chosen  practically  as  a basis 
for  modification  to  Eq.  (22).  The  early  stage  of  droplet  growth 

after  nucleation  takes  place  under  free-molecule-fiow  conditions  due  to  the 
small  cluster  size.  The  droplet  growth  rate  then  is  given  by  Eq.  (47a). 
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The  set  of  equations,  Eqs.  (67)  to  (71),  subject  to  stagnation 
conditions  completes  the  gas dynamic  description  of  the  flow  with  condensation. 
The  integration  of  these  equations  along  a atfreaioline^i,**  i.e.,  the  nozzle 
axis,  gives  the  flow  variables,  p,  u,  p,  T and  g (Fig.  9).  Provided  that  the 
pressure  ratio  across  the  nozzle  is  enough  to  produce  sonic  flow  at  the  throat 
and  that  thermal  choking  due  to  heat  release  of  condensation  does  not  occur, 
the  integration  is  rather  straightforward.  When  one  of  the  variables  such  as 
the  static  pressure  p(x)  is  measured  and  known,  one  of  the  remaining  variables, 
e.g.,  the  mass  fraction  of  condensate  g,  can  be  predicted  by  these  equations. 

To  test  the  nucleat ion- rate  equation  or  the  droplet- growth  law,  the  predicted 
variable  must  be  compared  with  the  measured  one.  Alternatively,  the  critical 
supersaturation , which  is  estimated  from  the  behaviour  of  the  condensate  mass, 
can  be  compared  with  the  value  calculated  using  the  nucleation-rate  equations. 

Employing  this  method  with  measurement  of  the  static  pressures  for 
vapours  of  HgO,  CgH^,  CHCI3,  CCI3F  and  C2H5OH  expanding  in  a supersonic  air 
nozzle  (wq  = 0.04  ~ 0.2),  Jaeger  et  al  (19°9)  and  Dawson  et  al  (1969)  concluded 
that  the  nucleation  rates  in  vapour  of  regular- shaped,  non-polar  molecules 
(NH3,  CgHg,  CHCI3,  CCI3F)  were  in  quantitative  agreement  with  the  Lo the- Pound 
theory,  Eq.  (22),  while  fluids  having  hydrogen  bonding  or  rod- shaped  molecules 
(H2O,  C2H5OH)  tended  to  follow  the  classical  theory,  Eq.  (18) . The  condensation 
coefficients  were  estimated  from  the  shape  of  the  pressure  profile  as  being 
between  1 and  0.1  for  the  former  fluids  ( NH3, C5H6 , CHCI3,  CCI3F)  and  between 
0.1  and  0.01  for  the  latter  (H2O,  C2H5OH) . The  orientation  of  molecules  at  the 
cluster  surface  is  considered  to  cause  these  feat  Tires  (see  Section  2.2). 

On  the  other  hand,  Wegener  et  al  (1972)  measured  the  static-pressure 
distribution  on  the  centreline  of  the  supersonic  nozzle  flow  of  ethanol-air 
(or  nitrogen)  mixture  with  a hypodermic  tube,  and  examined  a correction  factor 
to  the  classical  nucleation  theory  by  introducing  the  factor  r in  the  expression 
of  the  nucleation  rate 

1 - r • ici  (72) 

where  I c«  is  given  by  Eq.  (18).  Since  the  critical  radii  of  the  clusters  are 
appreciably  smaller  than  those  observed  in  cloud  chambers  owing  to  higher 
cooling  rates,  the  effect  of  droplet  radius  on  the  surface  tension  must  be 
taken  into  account  in  Ic£.  The  factor  r can  be  written  in  the  Lothe-Pound 
egression,  Eq.  (22), 
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By  calculating  the  steady- state  rate  of  nucleation,  the  replacement  factor 
qrep  and  the  effective  surface  tension  were  determined  so  as  to  match  the 
observed  condensation.  R value,  fpn  r can  be  chosen  for  which  the  rapid  increase 
in  calculated  nucleation  rate  coincides  with  the  observed  onset  of  condensation 
(Fig.  9).  When  the  surface  tension  takes  the  value  for  a flat  surface  vm  the 
experiment  yields  q^ep  "*  10^*  and  r **  105  for  Uq  - 0.05.  If  smaller  values  of 
surface  tension  are  chosen  for  smaller  droplets,  Qrep  becomes  larger  than  10l^ 
and  r approaches  unity;  for  o/er*  *0.8,  q^  **  lO®  and  r ■»  10.  A combination 
of  the  nucleation  rate  with  these  values  or  qrep  and  a and  the  droplet-growth 
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law  with  a suitable  condensation  coefficient  (ac  = 0.7)  permits  the  calcula- 
tion of  the  condensate  mass  fraction  to  fit  the  experimental  results  (Fig.  9) • 
Although  the  Lothe-Pound  factor  i^p  is  not  directly  applicable  to  liquid,  its 
values  (fLP  «■  1<A7)  axe  appreciably  above  those  from  experiments. 

Wu  et  al  (1977)  studied  condensation  of  sulfur  hexafluoride  (SFg) 
carried  in  argon  and  obtained  the  values  of  10*  £ r £ lO^  for  0.37  ^ Uq  ^ 0.01. 

As  mentioned  later,  using  rarefaction  waves  in  a shock  tube  for  water  vapour 
condensation  in  carrier  gases  of  He,  Ar  and  air  and  canparing  the  results  with 
the  non- isothermal  classical  theory,  Eq.  (7*0 , Barschdorff  (1975)  obtained 
10®  S r £ 10°  for  1 ^ Uq  ^ 0.004.  Wegener  and  Wu  (1976,  1977)  ccopared  these 
experiments  and  showed  the  values  of  r for  various  vapours  (Table  2).  All 
data  exhibit  higher  nucleation  rates  of  several  orders  than  those  predicted 
by  the  classical  theory.  It  is  interesting  to  compare  these  results  of  nozzle 
flows  r » 1 with  those  of  cloud  chambers  r *»  1,  although  the  reasons  for 
such  large  variations  remain  to  be  solved.  The  factor  r must  be  considered 
as  a correction  factor  for  the  classical  hcmogeneous-nudeation  rate  to  account 
for  both  contributions  from  statistical-mechanic  and  thermodynamic  sources  and 
from  other  uncertainties.  The  latter  may  arise  from  the  possible  non-hamogeneous 
nucleation  caused  by  inpurities  and  flow  properties. 

In  pure-vapour  condensation,  some  of  the  vapour  molecules  colliding 
with  droplets  must  act  as  a carrier  of  the  thermal  energy  from  the  latent  heat 
condensation.  Thus,  the  sinplifying  condition  of  equal  temperatures  for  drop- 
lets and  vapour  ceases  to  be  valid,  and  the  classical  isothermal  theory  will 
predict  considerably  higher  condensation  rates  or  smaller  values  of  r.  Intro- 
ducing a correction  factor  for  the  non-isothermal  removal  of  latent  heat  from 
condensate  clusters  £,  Barschdorff  et  al  (1971,  1976)  and  Ludig  (1975)  modified 
Eq.  (72)  for  wet  steam, 


1 - 5 r icJ  OHO 

The  factor  £ as  derived  by  Feder  et  al  (1966) , is 

C Lli (75) 

Comparing  the  calculated  pressure  distribution  with  the  measured  pressure 
profile,  they  obtained  the  best-fit  values  for  the  factor  r for  a number  of 
steam  experiments  from  several  investigations.  The  result  showed  that  classi- 
cal nucleation  theory  with  a correction  factor  £,  a bulk- liquid  surface-tension 
and  a thermal- accommodation  coefficient  of  unity  predict  well  the  experimental 
results  fof  low  initial  stagnation  pressures  (<0.5  bar) . For  higher  pressures, 
the  values  of  r decrease  up  to  about  ~ 10-5  at  5 bar . Saltan ov  et  al  (1971) 
introduced  a correction  factor  (3  for  nucleation  work  into  the  exponent  as 
exp(-0tf}*/ttO  and  compared  the  predicted  and  experimental  results  of  steam. 

The  factor  increases  with  the  stagnation  pressure,  taking  values  between  2.5 
and  5 in  the  range  of  1 < p0  < 32  bar.  This  also  means  smaller  values  of  r 
at  higher  pressures.  Although  quantitatively  unexplained,  one  of  the  main 
reasons  for  this  discrepancy  can  be  attributed  to  the  use  of  a perfect- gas 
equation  of  state . A real-gas  equation  reduces  the  nucleation  rate  considerably 
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(Bier  et  al,  1977).  Real-gas  analyses  were  made  by  Ludig  (1975),  Puzyrewski 
(1975),  Puzyrewski  and  StudzinskL  (1977)  and  Bier  et  al  (1977). 

When  coagulation  of  droplets  is  neglected,  the  number  of  droplets 
can  be  compared  directly  with  the  measured  value.  The  former  is  given  by 


n(x)  = p(x)  J*  4{  (76) 

o 

The  latter  can  be  obtained  by  means  of  light- scattering  methods  (Stein  and 
Wegener,  1967;  Stein,  1969;  Clumpner,  1971;  Stein  and  Moses,  1972;  and  Conard, 
1977).  By  this  technique,  Stein  and  Wegener  (1967)  found  agreement  for  water 
vapour  with  the  classical  nucleation-rate  theory  using  properties  of  ice  and 
a condensation  coefficient  of  unity.  Using  an  ethyl  alcohol -vapour  air  mixture, 
Clumpner  (1971)  studied  the  number  and  size  of  droplets  in  the  liquid  state 
leaving  a supersonic  nozzle  exit  and  showed  a narrow  distribution  of  droplet 
radius  about  a mean  size.  This  implies  that  most  of  the  droplets  formed  in  a 
narrow  range  of  locations  close  to  the  peak  value  of  the  nucleation  rate. 

In  actual  nozzle  flows,  two-dimensional  flow  conditions  between  the 
nozzle  axis  and  the  wall  exist  which  yield  different  cooling  rates.  Such  super- 
sonic-nozzle flows  with  condensation  were  studied  by  Davydov  (1972),  Tkalenko 
(1972)  and  Bratos  and  Jaescke  (1974)  using  the  method  of  characteristics.  The 
results  show  that  condensation  occurs  at  the  wall  further  upstream  than  at  the 
axis.  The  pressure,  temperature  and  density  deviate  more  strongly  from  the 
isentrcpic  values  along  the  axis  than  the  wall.  Also  the  distributions  of 
droplet  size  are  significantly  different  along  the  axis  and  the  wall.  These 
aspects  well  describe  the  experimental  results,  although  the  effects  of  two- 
phase  turbulent  motions  of  gas  and  droplet  and  of  the  wall  boundary  layer  were 
not  solved  in  the  analysis. 

Oscillating  unsteady  flows  are  often  encountered  in  supersonic-nozzle 
flows  with  condensation.  The  frequency  range  of  the  pressure  fluctuations 
extends  frcm  about  103  to  1C)4  Hz  and  the  amplitude  up  to  ± 5$  of  the  mean  value. 
These  phenomena  were  studied  theoretically  and  experimentally  by  Zierep  and  Lin 
(1968),  Barschdorf  and  Fillipov  (1970),  Barschdorf  ( 19^7,  1971),  Wegener  and 
Cagliostro  (1973)  and  Saltanov  (1974)  • The  cause  is  a feedback  system  caused 
by  the  movement  of  shock  waves  generated  by  the  condensation.  Such  aspects 
can  also  be  well  predicted  by  the  method  of  characteristics. 

From  these  results  it  may  be  said  that  the  features  of  nozzle  flows 
with  homogeneous  nucleation  and  condensation  can  be  predicted  by  one-  or  two- 
dimensional,  inviscid  gasdynamic  theory  coupled  with  a macroscopic  nucleation- 
rate  theory  and  a droplet  growth  law,  provided  that  an  empirical  adjustment 
factor  is  found  for  the  nucleation  theory.  In  order  to  physically  understand 
and  predict  the  adjustment  factor,  further  theoretical  and  experimental  studies 
of  the  nucleation  process  are  required  on  the  basis  of  molecular  dynamics  for 
clusterization . 

Applying  the  cluster  kinetic  theory  of  Buckle  (1969)  to  condensation 
of  water  vapour  in  air,  Pouring  (1977)  calculated  the  course  of  homogeneous 
condensation  in  a supersonic  nozzle  (see  Section  2.4) . It  was  shown  that  if 
a valid  initial-equilibrium  condition  is  well-established  with  the  assumed 


molecular  parameters , the  calculated  nonequilibrium  expansion  closely  matches 
the  experimental  pressure  distribution.  Although  the  estimation  of  the  mole- 
cular parameters  contains  sene  uncertainties,  this  approach  provides  a kinetic 
aspect  of  the  nucleation  and  condensation  process. 

Nozzle  flows  with  heterogeneous  condensation  are  often  observed  in 
chemical- engineering  processes  such  as  condensers,  filters,  and  cooling  towers. 
However,  owing  to  the  complications  introduced  by  heterogeneous  condensation, 
systematic  fundamental  studies  of  these  flows  as  yet  have  not  been  done.  Buckle 
and  Pouring  (1965)  studied  the  effect  of  large  numbers  of  smoke  particles  of 
silver  iodide  and  other  substances  (0.1  ~ 1 p diameters)  artificially  intro- 
duced in  a supersonic  nozzle  and  observed  no  difference  in  the  wall  pressure 
distributions  with  and  without  seeding.  Wang  (1976)  investigated  the  effect 
of  foreign  nuclei  on  nitrogen  condensation  with  a simple  model  that  they 
provide  only  additional  condensation  surfaces.  The  result  could  not  success- 
fully explain  the  experiment  showing  low  critical  super saturations . Saltanov 
and  Tsirihin  (1977)  studied  theoretically  heterogeneous  nucleation  of  nozzle 
flow  and  liquid- type  nuclei.  It  is  equivalent  to  homogeneous  nucleation  which 
has  droplets  of  a specified  size  and  concentration.  Condensation  of  molecules 
occurs  simultaneously  on  these  nuclei  and  those  nucleated  homogeneously.  It 
was  shown  that  the  addition  of  initial  droplets  by  more  than  0.1 $ of  the 
concentration  of  nuclei  nucleated  homogeneously  did  not  cause  the  pressure  to 
deviate  from  their  isentrqpic  values. 

For  heterogeneous  nucleation  and  condensation,  the  size-distribution 
of  nucleus  particles,  the  geometrical,  physical  and  chemical  characteristics 
of  the  nucleus  surface  and  the  electrical  interaction  between  the  nucleus 
particles  anu  the  gas  mixture,  take  on  Important  roles.  These  effects  must  be 
studied  to  ascertain  the  features  of  nozzle  flows.  At  high  droplet  population 
densities,  their  coagulation  due  to  Brownian-motion  or  turbulent- fluctuations 
will  be  appreciable.  Wegener  et  al  (1972)  estimated  that  (within  other  uncer- 
tainties) coagulation  could  have  little  effect  on  the  condensation  process  in 
a supersonic-nozzle  flow  of  an  ethanol-air  mixture.  However,  since  the  coagu- 
lation process  itself  has  unknown  complicating  factors,  more  quantitative 
studies  of  this  effect  are  required. 

Relative  motions  and  relative  differences  of  temperature  between 
the  condensate  droplets  : and  * the  gas  mixture  sire  usually  assumed 
negligible  (Wegener,  1969) . Under  seme  conditions  with  rapid  changes  in 
axial  velocity  and  temperature,  nonequilibrium  flows  result,  where  the  thermal 
and  kinematic  parameters  of  the  gas  and  the  droplets  differ.  This  leads  to 
internal  mass , momentum  and  energy  fluxes  between  the  constituents . The 
coupling  of  the  two  flow  systems  of  vapour  and  droplets  must  be  explicitly 
taken  into  account.  Such  studies  have  been  carried  out  by  Konorski  (1976)  for 
nozzle  flows  with  slight  deviations  from  the  equilibrium  state.  Further 
extension  of  these  studies  to  nucleation  and  condensation  processes  including 
particle  coalescence  would  extend  our  knowledge  of  flow  nucleation  and  conden- 
sation. 

4.3  Rarefaction  and  Expansion-Wave  Flows 

Nonstationeury  rarefaction-wave  flows  in  the  driver  section  of  a 
shock  tube  or  in  a Ludwieg  tube  or  steeidy  expansion-wave  flows  around  a sharp 
corner  or  in  nozzles  can  result  in  vapour  condensation,  end  offer  additional 
techniques  for  the  study  of  nucleation  and  condensation.  Earlier  work  on  these 
subjects  up  to  1974  v&a  reviewed  by  Wegener  (1975). 
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Steady-expansion  flow  around  a sharp  comer  (a  Prandtl-Meyer  expansion, 
Fig.  10)  provides  cooling  rates  which  decrease  strongly  away  frcm  the  comer. 

The  cooling  rate  varies  from  infinity  at  the  comer  to  low  values  at  large 
distances  from  the  corner.  Hence,  for  flows  of  condensing  vapour , the  features 
of  nucleation  and  condensation  depend  appreciably  on  the  radial  distance  from 
the  corner.  Since  higher  cooling  rates  result  in  higher  supercooling  at  the 
onset  of  condensation,  the  onset  supercooling  decreases  with  distance  from  the 
corner.  These  phenomena  were  studied  theoretically  and  experimentally  by 
Steffen  (1967) , Tkalenko  (1970),  Smith  (1971),  Kurshakov  et  al  (1971),  and 
Frank  (1977) • Smith  (1971)  studied  experimentally  the  expansion  of  moist  air 
around  a 40*  sharp  corner  which  was  part  of  a lower  block  of  a two-dimensional 
nozzle.  A Mach-Zehnder  interferometer  was  used  to  measure  the  density  distri- 
bution. The  result  showed  that  the  condensation  zone  was  a strong  function  of 
both  expansion  angle  e and  radial  distance  r.  from  the  corner,  as  expected. 

The  condensation  front  is  curved  (concave)  to  the  oncoming  flow  and  becomes 
parallel  to  the  equilibrium  characteristic  behind  the  head  of  the  expansion  fan 
with  increasing  distance  frcm  the  corner.  The  super cooling,  Ts  - Tc,  where  Ts 
is  the  saturation  temperature  and  Tc  the  temperature  at  the  onset  of  condensa- 
tion, can  be  expressed  as  a function  only  of  the  local  cooling  rate  (-  dT/dt): 

Tg  - Tc  - 40.8  (t:  nsec)  (77) 

for  the  range  of  1.7  > (-  dT/dt)  > 0-3  K/psec.  Frank  (1977)  also  investigated 
such  flow  fields  of  higher  humidities  (uq  > 0.01)  and  observed  oblique  shock 
waves  within  the  expansion  fan.  At  higher  humidities,  two  clearly  separated 
condensation  zones  were  observed. 

Barrand  and  Rieutord  (1973)  studied  experimentally  the  nonstationary 
rarefaction-wave  flow  of  water  vapour-argon  or  helium  mixtures  in  a long  duct 
(8  cm  dia  x 201  cm  long)  The  pressure  history  was  measured  at  three  points 
along  the  duct  to  detect  the  onset  of  condensation.  With  these  results  and 
other  thermodynamic  variables  obtained  fran  measurements  of  dry-gas  expansions, 
they  calculated  a relation  between  cooling  rate  and  supercooling  at  the  onset 
of  condensation  which  was  approximately  an  extension  of  Smith's  results  (Fig. 

11) : 

✓ v0.12 

Ts  - Tc  = 0 ~ 2)  x 40.8  (-  § J (t:  psec)  (78) 

for  0.1  > (-  dT/dt)  >0.02  K/psec. 

In  addition  to  cloud- chanber  and  nozzle- flow  methods,  the  well- 
established  shock-tube  technique  provides  one  more  possibility  for  investiga- 
ting experimentally  and  analytically  nucleation  and  condensation  processes. 

A shock  tube  in  its  simplest  form  consists  of  a straight  constant-area  duct 
which  is  divided  by  a removable  diaphragm  into  a high-pressure  driver  section 
and  a low-pressure  driven  section.  By  bursting  the  diaphragm,  a shock  wave 
rapidly  develops  and  moves  into  the  driven  section,  and  an  i sen tropic  rare- 
faction wave  propagates  into  the  driver  section  (Fig.  12).  Since  the 
gas dynamic  behaviour  of  such  flows  is  now  well-understood  and  easily  controlled, 
rarefaction  waves  have  proved  most  useful  for  investigating  condensation.  Other 
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applications  of  shock  tubes  to  condensation  studies  were  reported  by  Rung  and 
Bauer  (1971),  Hamer  et  al  (1971),  and  Hamer  and  Hurle  (1972). 

Rarefaction  waves  will  suddenly  expand  isentropically  mixtures  of 
condensing  vapour  and  noncondensible  gas  initially  at  rest  in  the  driver 
section.  The  mixtures  are  cooled  up  to  or  below  the  super saturation  state. 
Further  expansion  leads  to  appreciable  condensation  of  the  vapour,  which 
diminishes  the  supersaturation  state.  As  the  rarefaction-wave  head  moves  into 
the  mixture  at  rest  in  the  driver  with  the  state  sound  speed,  the  supersatura- 
tion along  a particle  path  reaches  a maximum  and  substantial  cluster  formation 
occurs.  The  onset  of  condensation  can  be  observed  by  a sudden  pressure  rise 
with  respect  to  the  isentrqpic  expansion  and/or  abrupt  increase  in  the  light 
scattering.  Since  in  a flow  with  a higher  cooling  rate  the  onset  of  condensa- 
tion occurs  at  a higher  super cooling,  the  first  onset  in  a shock  tube  takes 
place  close  to  the  tail  of  the  rarefaction  fan  at  a point  close  to  the  diaphragm 
where  the  cooling  rate  is  high.  Far  from  the  diaphragm  the  onset  moves  toward 
the  head  of  the  rarefaction  fan.  Thus,  the  condensation-wave  front  has  a 
curved  path,  as  shown  in  the  (x,  t) -diagram  or  distance- time  plane  (Fig.  13) • 

The  heat  released  by  condensation  causes  an  appreciable  change  in  the  flow 
variables.  Rarefaction  waves  within  the  driver  section  induce  a subsonic  flow 
behind  them.  The  heat  release  changes  the  gasdynamic  state  along  character- 
istics from  the  usual  isentropic  values  downstream  of  the  condensation  zone, 
as  shown  by  the  bump  in  the  (p,  x)  and  (T,  x)  sketches  at  the  "onset"  (Fig.  13)- 
Inside  the  condensation  zone  a two-phase  flow  of  droplets  and  gas  takes  place, 
and  under  appropriate  conditions,  shock  waves  are  produced  (Sislian  and  Glass, 
1976) . 


Actual  rarefaction  waves  can  be  viewed  effectively  as  a centred-wave 
system  with  a virtual  origin  not  far  from  the  actual  diaphragm  location.  The 
virtual  origin  is  determined  by  effects  of  fluid  inertia  for  the  initial 
acceleration,  characteristics  of  diaphragm  rupture  and  three-dimensional 
effects  due  to  diaphragm  curvature  and  associated  shock- wave  formation.  Thus, 
close  to  the  diaphragm,  ideal  simple  centred  waves  are  not  realized.  At  large 
distances  from  the  diaphragm,  the  ideality  also  deteriorates  owing  to  the  growth 
of  wall  boundary  layers.  Once  the  virtual  origin  is  determined  from  measurements 
of  static  pressure  the  rarefaction  fan  in  the  intermediate  distance  range  can 
be  predicted  frcm  one- dimensional  inviscid  gasdynamics  (Hall  et  al,  197*0  • 

The  shock-tube  technique  has  the  advantages  of  providing  Well- controlled 
experimental  conditions  which  have  the  combined  features  of  the  cloud  chamber 
and  the  supersonic  nozzle.  Like  the  cloud  chamber,  only  small  amounts  of 
condensing  substances  are  required,  hence  expensive  or  toxic  species  can  be 
studied.  It  is  also  feasible  to  control  vapour,  contamination  by  pimping  to 
high  vacuum  prior  to  the  experiments . As  for  supersonic-nozzle  flow,  shock- 
tube  operation  is  well  controlled  and  understood  from  one- dimensional  inviscid 
gasdynamics,  provided  the  initial  pressures  (Reynolds  numbers)  are  not  too  lew. 

The  flow  history  along  a particle  path  differs  depending  on  its  distance  from 
the  diaphragm.  At  different  locations  along  the  tube,  a wide-range  of  cooling 
rates  is  obtainable  frcm  those  of  cloud  chambers  (10“2  K/|is)  to  those  of  super- 
sonic nozzles  (1  K/ns)  (Table  1) . Usually,  the  cooling  rate  is  lower  than  in 
nozzle  flows.  However,  contrary  to  cloud  chamber  and  supersonic  nozzle,  the 
flow  in  shock  tubes  is  nonstationary.  During  rapid  expansion  (near  a corner  or 
a diaphragm)  nucleation  takes  place  in  a nonsteady  process  and  the  steady-state 
nucleation  theory  may  not  apply.  However,  for  the  range  of  usual  cooling  rates 
mentioned  above,  quasi- steady- state  nucleation  may  be  established  in  the  available 
flow  times. 
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Recently,  Kawada  and  Mori  (1973)  successfully  applied  the  rarefaction 
wave  to  study  the  nucleation  and  condensation  rates  of  several  refrigerant 
vapours  by  measuring  the  static  pressure,  density  (Mach-Zehnder  interferometer) 
and  light  transmission.  The  onset  of  condensation  was  defined  as  the  point 
where  the  transmitted  light  was  attenuated  by  %.  The  transmission  efficiency 
of  light  is  related  to  the  size-distribution  of  droplets  which  can  be  predicted 
from  nucleation  and  condensation  rates.  Thus,  by  measuring  the  onset  time  at 
a fixed  location  in  the  tube,  the  nucleation  theory  can  be  tested  by  experiments. 
They  showed  that  the  classical  theory  predicts  well  the  experimental  results 
with  values  of  surface  tension,  (r/<T„  = O.38  for  R-21,  O.69  for  R-4o,  0.6l  for 
R-113  + argon,  and  £ 1 for  CH3OH  + argon 

* 

Using  rarefaction  waves,  Barschdorff  (1975)  also  investigated  the 
nucleation  rate  of  water  vapour  in  different  inert  gases  (air,  helium  and 
argon)  over  a wide  range  of  mixing  ratios,  and  studied  the  carrier-gas  effects 
on  homogeneous  nucleation  of  water  vapour  The  shock  tube  was  evacuated  to 
10-6  torr  using  a diffusion  pump.  Distilled,  purified  water  vapour  and  a filt- 
ered, ultra-high-purity  carrier  gas  (argon,  helium  or  dry  air)  filled  the  driving 
chamber.  The  rarefaction  process  was  studied  by  measuring  the  time  histories 
of  static  pressure  and  light  scattering.  As  in  supersonic-nozzle  flows,  the 
mass  fraction  of  condensate  g was  computed  by  integrating  Eq.  (55)  along  a 
particle  path  with  a macroscopic  nucleation  theory,  Eq.  (74) • When  the  onset 
of  condensation  is  of  primary  interest,  the  gasdynamic  state  of  flow  can  be 
greatly  simplified  by  adjusting  the  initial  diaphragm  pressure  ratio  so  as  to 
make  the  onset  coincide  with  the  tail  of  the  rarefaction  wave  at  the  observation 
point.  Then,  there  is  no  influence  of  prior  condensation  in  the  rarefaction 
fan  and  an  isen tropic  state  is  preserved  Tip  to  the  onset  (Wu,  1977) . The  value 
of  r in  Eq.  (74)  was  chosen  for  which  the  strong  increase  in  nucleation  rate 
matched  well  with  the  observed  onset  of  condensation.  The  results  show  a 
variation  of  1C)6  £ p £ 1 between  the  limits  of  low  vapour  concentration 
(wo  = 0.003)  to  pure  steam  (wq  = 1) . No  systematic  deviations  due  to  the  type 
of  carrier  gas  were  observed.  The  calculation  of  r was  performed  by  using 
the  surface  tension  corresponding  to  the  bulk  value  cr,,,  and  the  condensation 
coefficient  of  unity  in  the  droplet  growth.  Although  this  assumption  may  explain 
partially  the  large  deviation  of  the  r values,  further  theoretical  work  is  neces- 
sary for  its  understanding  on  a basis  of  two- component,  two-phase  statistical 
mechanics  or  kinetics. 

Lee  (1977,  1978)  employed  the  same  method  for  evaluating  his  experi- 
ments on  condensation  of  water,  heavy  water  (D2O) , benzene  (C^Hg) , carbon 
tetrachloride  (CClij'),  and  freon  11  (OCI3F)  in  argon.  No  difference  was  obser- 
ved in  the  condensation  of  H2O  and  D2O . The  values,  of  r calculated  were  103 
~ 105  for  H2O,  10l  ~ 103  for  D20  and  C6H£,  102  ~ 1CT  for  CClj.,  and  109  ~ 1012 

. for  freon  11. 

Kalra  (1975)  and  Glass  et  al  (1977)  investigated  the  condensation  of 
water  vapour  in  a nitrogen  carrier  gas.  The  time  variations  of  density  gradient 
and  static  pressure  as  well  as  the  onset  of  condensation  were  monitored  at  two 
fixed  locations  in  the  driver  section  by  using  a laser  Fabry-Perot  interfero- 
meter, a differential  interferometer  and  piezotron  transducers.  The  effects 
of  the  cooling  rate  on  the  super saturation  at  the  onset  of  condensation  provided 
an  empirical  relations 
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= 1.73  x 10 


(79) 


for  0.09  > (-  dT/dt)  > 0.04  jisec.  This  result  coincides  approximately  with 
those  of  Smith  (Eq.  77) , and  Barrand  and  Rieutord  (Eq.  78),  as  seen  in  Fig.  11. 

Several  analyses  of  rarefaction- wave  flow  of  a mixture  of  water  vapour 
and  nitrogen  in  a shock  tube  were  done  by  Sislian  (1975) , Sislian  and  Glass 
(1976),  Kotake  and  Glass  (1976,  1977),  Glass,  Kalra  and  Sislian  (1977)  and  Glass 
and  Sislian  (1977)*  The  method  of  characteristics  and  the  Lax  method  were  used 
along  with  a macroscopic  nucleation  theory.  To  construct  a mathematical  formula- 
tion of  one-dimensional  nonstationary  flow  in  a shock  tube,  the  same  assumptions 
vfere  employed  as  for  nozzle  flows:  an  inviscid,  nonconducting  perfect  gas 
mixture,  no-slip  velocity  condition,  identical  temperatures  of  condensate  and 
gas  mixture,  and  a uniform  distribution  of  the  condensate  across  the  tube 
section.  The  governing  equations  are: 


dh  _ 1 dp 

dE  p dK 


(80) 


with  t - time,  and  x - distance  from  the  diaphragm  along  the  shock  tube.  The 
equation  of  state  is  given  by 


and  the  specific  enthalpy  is 


h = Cpo1  ’ 

where 

°Po  ‘ (1  ‘ “o  ■ «p)op*  * “o  V 4 ®P  5* 


The  production  of  the  condensate  mass  along  a particle  path  can  be  expressed 
by  Eq.  (56)  or  Eq.  (63)  in  which  the  nucleation  rate  and  the  droplet-growth 
rate  are  given  by  Eqs.  (18)  or  (39)  and  (47),  respectively. 

The  above  set  of  equations  were  solved  numerically  using  a second 
order  characteristic  method  subject  to  the  initial  conditions.  Sislian  (1975) 
and  Sislian  and  Glass  (1976)  calculated  the  case  of  homogeneous  condensation 
for  water  vapour-nitrogen  mixture  of  uq  = 0.01776  and  T0  = 295.3  K (Figs.  13, 
l4  and  15) . As  shown  in  Fig.  13,  the  characteristics  emanating  from  the 
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origin  remain  straight  as  long  as  the  effect  of  homogeneous  condensation  on 
the  flow  variable  is  negligibly  small.  After  penetrating  the  condensation 
zone  the  characteristics  begin  to  curve  appreciably  towards  the  head  of  the 
rarefaction  wave  due  to  the  heat  release  from  condensation.  Computations 
showed  that,  for  the  case  considered,  the  condensation  wave  is  followed  by 
a shock  wave  generated  by  the  intersection  of  characteristics  of  the  same 
family.  The  width  of  the  condensation  wave  (zone)  increases  with  increasing 
distance  from  the  origin  due  to  the  interplay  of  two  varying  opposing  effects 
of  decreasing  pressure  due  to  expansion  and  increasing  pressure  due  to  the 
heat  release  from  condensation.  Close  to  the  origin  where  the  expansion  rate 
is  very  high,  appreciable  condensation  takes  some  time  to  occur.  However,  in 
this  region,  the  steady-state  nucleation  theory  may  not  be  applicable  in  order 
to  obtain  detailed  flow  information. 

The  presence  of  shock  waves  in  the  flow  field  complicates  the  numerical 
solution  of  the  problem.  In  order  to  handle  such  discontinuities,  Lax's  method 
of  Implicit  artificial  viscosity  is  used.  This  technique  requires  that  the 
governing  equations  be  expressed  in  Lagrangian  form: 


(81) 


with  5 - Lagrangian  parameter  which  identifies  a fluid  element,  taken  as 
the  spatial  coordinate  at  the  time  origin  t = 0;  - density  of  the  fluid  element 

5 at  t = 0.  The  agreement  between  both  calculations  was  satisfactory  except  in 
the  region  close  to  the  origin.  The  analyses  were  conpared  with  the  experiments 
of  Karla  (1975)  and  Glass  et  al  (1977)  (Fig.  15) • The  pressure  history  shows  a 
significant  discrepancy  between  analysis  and  experiment . If  the  value  of  the 
macroscopic  surface  tension  (era,  = 82  dyne/cm)  is  adjusted  to  a lower  value 
(cr  = 68  dyne/cm),  = O.83,  the  analysis  can  be  made  to  match  the  experi- 

mental data  fairly  well  except  for  the  value  of  supersaturation  at  the  onset 
of  condensation. 


Kotake  and  Glass  (1976,  1977)  studied  an  alternative  heterogeneous 
condensation  analysis.  Fine  condensation  nuclei  of  solid  particles  (Aitken 
nuclei)  with  radii  below  0.1  pm  were  assumed  to  exist  after  usual  filtrations 
(Mason,  1971,  Boucher,  1969) . A size  distribution  of  nuclei  was  assumed 

A^ar-Kw1'5}  <« 

where  R0  ■ (1.5  ~ 5)  x 10"6  an  (Fig.  1$).  All  particles  were  assumed :to 
have  a smooth,  spherical  shape  of  effective  radius  without  any  coagulation  or 
separation  which  were  chemically  inert  and  insoluble  in  water  vapour.  Rare- 
faction wave  flows  with  heterogeneous  condensation  were  numerically  calculated 
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using  macroscopic  heterogeneous-nucleation  theory,  Eq.  (39) > and  condensate- 
mass  equation,  Eq.  (63),  coupled  with  the  flow  equations.  Owing  to  the 
complexity  and  variety  of  collisions  between  monomers  and  the  solid  surface 
of  nuclei,  it  is  difficult  to  estimate  the  values  of  the  vibrational  frequen- 
cies of  atoms,  the  desorption  free-energy  of  monomers  and  the  contact  angle. 
Fortunately,  the  first  two  factors  are  not  significant  compared  with  the 
contact  angle.  The  computations  showed  that  the  features  of  heterogeneous 
condensation  are  dominated  largely  by  the  contact  angle. 

The  contact  angle  directly  affects  the  activation  energy  .The  ydluft  f=l8o* 
implies  the  same  value  of  activation  energy  as  homogeneous  nucleatian,  whereas 
8 = 90®  yields  a value  reduced  by  about  50%.  Less  activation  energy  results 
in  a higher  nucleation  rate  with  lower  super saturation,  lower  supercooling,  a 
higher  propagation  velocity  of  the  condensation  wave  and  a narrower  condensation 
zone.  The  calculated  results  show  these  patterns  by  comparison  with  the  homo- 
geneous case.  By  choosing  a suitable  value  for  the  contact  angle  (e  »90*)>  the 
experimental  data  could  be  fitted  fairly  well,  again  except  for  the  super  satura- 
tion at  the  onset  of  condensation  (Fig.  15). 


By  using  macroscopic  models  of  homogeneous  and  heterogeneous  nucleation 
with  a suitable  value  of  the  surface  tension  and  the  average  contact  angle, 
respectively,  condensation  in  rarefaction  waves  in  a shock  tube  can  be  predicted. 
Although  the  experiment al  results  can  be  explained  equally  well  by  either  homo- 
geneous or  heterogeneous-nucleation  theory,  many  factors  remain  unsolved  such 
as  coagulation,  solubility  and  electrical  charges.  Nevertheless,  the  procedure 
of  choosing  values  of  surface  tension  or  contact  angle  to  fit  experimental 
data  is  not  a satisfactory  one.  At  best  they  are  a guide  for  decisive  quantita- 
tive measurements  of  these  parameters  under  actual  conditions. 

The  vapours  treated  above  can  be  condensed  to  the  liquid  phase  only 
by  expansion.  However,  the  so-called  "retrograde"  fluids  can  change  their 
phase  from  vapour  to  liquid  by  "compression".  In  the  phase  change,  these 
fluids  store  the  latent  heat  in  the  condensate  state  within  a large  number  of 
internal  (vibrational)  molecular  degrees  of  freedom.  This  is  the  case  if  the 
heat  capacity  of  a fluid  at  a reference  state  is  sufficiently  large.  All  the 
fluids  in  the  paraffin  family,  say  hydrocarbon  series  above  the  complexity 
level  of  propane,  n-hexane  or  n-octane,  show  this  retrograde  behaviour  and 
under  some  conditions  condensation  from  the  vapour  phase  to  the  liquid  phase 
will  be  produced  by  a shock  wave. 

Thompson  and  Sullivan  (1975)  studied  the  possibility  of  such  "lique- 
faction" shock  waves.  They  showed  that  for  fluids  which  have  an  ideal  gas 
specific  heat  Oy  > 2k  R at  the  critical  tenperature,  "complete  condensation 
or  liquefaction  shocks"  can  be  found  in  relatively  simple  experiments.  A 
"complete  condensation  shock"  means  that  the  upstream  state  is  "fully  dry" 
and  the  downstream  state  is  "fully  wet".  Dettleff,  Thompson  and  Meier  (1976) 
and  Dettleff  et  al  (1976)  carried  out  experimental  studies  of  condensation 
induced  by  liquefaction  shocks  using  a shock  tube  with  various  fluorinated 
hydrocarbons  as  the  test  fluid.  Measurements  of  pressure,  shock  velocity  and 
temperature  showed  good  agreement  with  the  Rankine-Hugoniot  relations,  and  a 
clear  liquid-phase  condensate  was  photographically  observed  behind  the  shock 
front.  Although  further  studies  are  required  for  a quantitative  description 
of  the  flow  field,  this  new  aspect  of  condensation  will  provide  additional 
approaches  to  the  development  of  nucleation  and  condensation  theories. 


4.4  Free-Molecular  Jets 

Most  experiments  on  nucleation  and  condensation  of  vapour  were 
limited  to  examining  the  supersaturation  at  which  appreciable  condensation 
occurs.  Only  information  integrated  over  many  complex  intermediate  steps  is 
available.  Consequently,  the  details  of  the  collision  kinetics  are  hidden. 

A quantitative  interpretation  of  each  step,  from  monomer  or  dimer  collisions 
to  droplet  formation,  is  difficult  to  obtain  from  macroscopic  measurements 
such  as  droplet  sine,  pressure,  density  or  temperature.  To  provide  all  the 
desired  information  on  nucleation,  an  ideal  experiment  would  make  it  possible 
to  observe  each  stepwise  event.  Molecular-beam  saapling  of  a free  jet  seems 
capable  of  providing  such  information.  For  this  reason,  it  has  been  used 
intensively  during  the  past  decade  for  the  study  of  nucleation  kinetics  of 
early  clustering  stages. 


A free-jet  expansion  consists  of  allowing  a gas  in  colli si  on- 
dominated  continuum  flow  to  expand  through  an  orifice  or  nozzle  into  a region 
which  is  maintained  at  high  vacuum  (Fig.  17).  As  the  gas  expands  from  the 
source  into  a vacuum,  random  thermal  motion  is  converted  to  directed  motion 
and  the  gas  is  accelerated  to  sonic  flow  at  the  nozzle  throat.  Beyond  this, 
the  gas  near  the  centre  of  the  jet  expands  supersonically  as  a nearly-isentropic 
continuum  flow  until  it  becomes  rarefied  enough  to  become  a free-molecular  flow. 
Typical  pressures  drop  three  to  four  orders  of  magnitude  and  temperatures  fall 
to  a few  degrees  Kelvin.  The  rate  of  expansion  is  controlled  by  the  orifice  or 
nozzle  diameter,  since  the  expansion  history  of  the  gas  is  scaled  with  distance 
(in  diameters)  downstream  of  the  sonic  throat.  As  shown  schematically  in  Fig. 
l6,  the  core  of  this  low-density  jet  is  then  usually  skimmed  into  a molecular 
beam  through  a small  hole  at  the  tip  of  a conical  "skimmer"  located  downstream 
of  the  nozzle.  Usually,  continuum  flow  exists  upstream  of  the  skimmer  and 
free-molecular  flow  downstream.  Under  many  conditions,  a transition  from 
continuum  flow  to  free-molecular  flow  occurs  upstream  of  the  skimmer.  Mole- 
cules flying  near  the  system  centreline  are  then  passed  tnrough  a collimating 
orifice  into  a detection  chamber. 

During  the  initial  high-density  stage  of  the  expansion,  the  preci- 
pitous adiabatic  decrease  of  tenperature  relative  to  pressure  may  initiate 
homogeneous  nucleation  if  the  gas  becomes  supersaturated  and  the  growth  of 
molecular  clusters  can  occur  until  the  collision  rate  decreases  for  molecular 
flow  (Fig.  l8) . If  the  sampling  process  is  carried  out  without  materially 
disturbing  the  flow,  one  can  obtain  a so-called  condensed  molecular  beam, 
i.e.,  cluster  beams  in  arbitrary  stages  of  condensation  arising  from  upstream 
stagnation  conditions.  They  are  made  up  of  clusters,  i.e.,  microdrqplets  or 
microcrystals  where  atoms  or  molecules  are  held  together  by  Van  der  Walls 
forces.  Clusters  are  viewed  as  species  with  properties  which  are  intermediate 
between  those  of  homogeneous  monomolecular  vapour  and  those  of  bulk  condensate. 
Determination  of  the  beam  composition  can  be  conveniently  carried  out  by  gas- 
scattering methods,  light- scattering  methods,  electric  or  magnetic  mass 
spectrometers,  retarding  mass  spectrometers,  time-of -flight  mass  spectrometers 
and  electron-beam  diffraction  methods. 

Figure  l8  shows  the  increase  in  Mach  number  M with  distance  x/d  from 
the  throat  to  its  design  value,  the  rapid  drop  in  tenperature  T and  the  increase 
in  concentration  of  the  condensate  clusters  n2»  no,  ...  . The  kinetic  process 
of  clustering  produces  clusters  of  smaller  size  at  higher  rates  (n2  > nj  >..».) 
of  which  concentrations  are  considerably  lower  than  the  equilibrium  values 
corresponding  to  the  tenperature. 
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In  the  last  decade,  the  study  of  dimer  and  cluster  formation  in 
supersonic  free  jets  has  been  advanced  by  many  investigators,  such  as:  Milne 
and  Greene  (1967,  E) , Milne,  Vandegrift  and  Greene  (1970,  A),  Golomb,  Good 
and  Brown  (1970,  E,  A),  Gordon,  Lee  and  Herschbach  (1971,  E,  A),  Golomb  et  al 
(1972,.  E),  Hagena  and  Obert  (1972,  E),  Dittmers,  Schutze,  Fischer  and  Schttgerl 
(1972,  E),  Yealland,  Deckers,  Scott  and  Tuori  (1972,  E),  Stein  and  Armstrong 
(1973,  E),  Sinha,  Schultz  and  Zare  (1973,  E),  Dorfeld  and  Hudson  (1973a,  b,  A), 
Van  Deursen  and  Reuss  (1973,  E),  Good  and.  Golomb  (1974,  E,  A),  Hagena  and 
Wedel  (1974,  E),  Armstrong  and  Stein  (1974,  E) , Seary  and  Fenn  (1974,  E), 

Beylich  and  Robben  (1974,  E,  A),  Seary  (1975,  A),  Calo  (1975,  E,  A),  Beylich 
(1976,  E,  A),  Van  Deursen  and  Reuss  (1977,  E),  Van  Deursen  and  Reuss  (1977,  E) , 
Khuth  (1977,  A),  where  E and  A denote  experimental  and  analytical  studies, 
respectively. 

In  Table  3 the  experimental  studies  are  listed,  showing  the  investi- 
gator, species  of  cluster,  diameter  of  nozzle  or  orifice  (d),  range  of  source 
pressure  (pQ)  and  temperature  (T0),  type  of  detector,  and  the  main  results. 

For  gases  such  as  Ar,  Kr,  Ne,  Xe,  N2,  ®2>  N0,  00,  SO2,  CO2,  K,  Rb,  Cs,  Na  and 
H2O,  cluster  concentrations  and  cluster  sizes  have  been  measured  mass-spectro- 
metrically  as  a function  of  source  pressure,  temperature,  and  nozzle  diameter. 
Figure  19  shows  typical  total  and  individual  mass  beam  intensities  against 
source  pressure.  As  the  source  pressure  is  increased  with  keeping  the  source 
temperature  constant,  the  total-beam  intensity  increases  slowly  to  a plateau 
due  to  the  larger  number  of  collisions  resulting  from  the  longer  continuum 
e3q)ansion  and  the  initial  state  closer  to  saturation-  After  the  plateau,  it 
follows  a slow  decrease  and  then  a steep  increase.  The  decrease  after  plateau 
is  attributed  to  skimmer  interaction.  The  steep  increase  is  due  to  heavier 
clusterization,  being  associated  with  decreasing  monomer  intensity  and  increasing 
dimer  and  cluster  intensities.  After  reaching  a maximum,  the  dimer  and  cluster 
intensities  decline  and  indicate  condensation  or  coalescence  of  clusters  into 
droplets  which  are  not  detected  by  the  mass  spectrometer.  Thus,  the  peaks  of 
dimer  or  cluster  curves  can  be  related  to  the  condensation  process . The  position 
of  the  peaks  depends  on  the  expansion  history,  that  is,  the  source  pressure  po» 
temperature  To  and  the  nozzle  or  orifice  diameter  (d)  in  the  form  d^PoT^11  = 
constant . 

For  correlation  of  the  onset  of  cluster  condensation,  the  following 
relations  have  been  obtained  experimentally  or  analytically  by  several  investi- 
gators: 


d°*55p  T-2*2  = const.  for  Ar,  Golonto  et  al  (1970) 
o o 


*0-6pot;2-2 


dO-W2-4 

o o 


= const.  for  Ar,  Kr,  COg,  SOg,  Ng,  CO  and  NO, 
Good  and  Golanfc  (1974) 

= const.  for  Ne,  Ar,  Kr,  Xe,  Ng  and  COg, 
Hagena  and  Obert  (1972) 


d°*55pcr"2*25  = const,  for  Ar,  Ne,  Hg,  Ng  and  0g, 

Van  Deursen  and  Reuss  (1977) 
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A more  detailed  knowledge  of  how  gas  properties  and  the  source  conditions 
including  nozzle  geometry  affect  cluster  formation  can  he  obtained  by  con- 
sidering similarity  conditions  of  the  phenomena.  Such  similarity  relation- 
ships have  been  developed  by  Hagena  (1969,  1974) , JIagena  and  Obert  (1972)  and 
Good  and  Golornb  (1974) . 


Fundamentally,  the  process  of  cluster  formation  may  be  solved  by 
using  the  appropriate  kinetic  equations  of  clusterization.  The  propriety  of 
the  kinetics  can  be  confirmed  by  comparing  the  results  with  experiments.  In 
order  to  arrive  at  meaningful  results,  a substantial  knowledge  is  required 
of  free- jet  expansion  from  continuum  flow  to  free-molecular  flow.  However, 
a complete  description  of  such  flow  fields  including  the  subsonic  region  has 
never  been  given  for  a small  orifice  or  nozzle.  It  is  not  possible  at  present 
to  give  a complete  history  of  the  states  to  which  the  clusters  may  be  exposed. 
During  the  very  first  stage  of  the  expansion,  the  state  of  the  system  changes 
slowly  enough  that  equilibrium  is  maintained.  As  the  expansion  proceeds,  the 
time  scale  for  change  in  the  clustering  environment  becomes  so  i that 
equilibrium  cannot  be  maintained.  Departure  from  equilibrium  occurs  as  soon 
as  the  flow  leaves  the  nozzle.  Experimentally  and  analytically,  it  is  shown 
4k&t  the  expansion  of  chemically  inert  gases  remains  isentrcpic  up  to  a so- 
called  "terminal  Mach  number"  which  can  be  calculated  from  the  gas  properties 
and  the  orifice  or  nozzle  size.  In  this  range,  using  numerical  methods,  the 
state  of  the  gas  in  a free  jet  can  be  computed  for  a given  orifice  or  nozzle 
and  the  stagnation  conditions.  Calculations  of  such  flows  have  been  carried 
out  for  various  specific-heat  ratios,  and  algebraic  fitting  relations  are 
available  which  provide  a close  approximation  to  the  numerical  results. 

Further  expansion  quickly  brings  about  a transition  to  free-molecular  flow. 

For  the  transition  process  from  continuum  to  free-molecular  flow,  a number  of 
studies  have  been  done  in  the  last  decade  by  using  methods  such  as  the  sudden- 
freeze  model,  simple  kinetic  theory,  BGK-model  and  Monte  Carlo  methods.  For 
large  Mach  numbers  (M  > 3)  and  large  nozzle  Reynolds  numbers  (Re  > 200) , the 
gas  remains  in  thermodynamic  equilibrium  up  to  the  terminal  Mach  number. 
Thereafter,  its  axial  translational-energy  distribution  and  the  parallel 
temperature  remain  "frozen"  at  the  terminal  Mach- number  value. 


Application  of  the  above  flow  features  to  the  clustering  process 
in  a free-jet  expansion  makes  it  feasible  to  analyze  the  kinetic  process  of 
cluster  formation  without  introducing  troublesome  and  detailed  flow  struc- 
tures. For  one-dimensional  isentrcpic  expansion  with  negligible  effects  from 
molecular  transports,  thermal  radiation  and  external  fields,  the  equations  of 
mass,  momentum  and  energy  give 


dx 

= 0 

dU  dn 

pudS+l£ 

= 0 

d / . 1 2 \ 

51  *+2U 

(83) 


where  p is  the  density  of  the  gas  and  cluster  mixture,  u the  axial  velocity, 
p the  pressure  and  h the  enthalpy . Although  the  cross-sectional  area  of  the 
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flow  Af  Is  an  important  variable  and  is  not  available  without  solving  the 
two-dimensional  problem,  the  result  for  a free- jet  expansion  of  chemically 
inert  gases  can  be  used  approximately.  The  equation  of  state  is 


p = pRT 


(84) 


where  R is  the  gas  constant  of  the  mixture.  The  rate  equation  of  cluster 
formation  can  be  written  as 


<*ni 

u = (ni»  •••»  “n*  P#  ^)  (85) 


where  ni  is  the  concentration  of  i-mer.  Generally,  the  enthalpy  is  expressed 
as  the  sum  of  i-mer  enthalpies 

h=y^('Vms+hrot+14lb+^1) 

where  htrans,  hrot  ^ Wib  denote  enthalpy  contribution  due  to  i-mer  trans- 
lational, rotational  and  vibrational  motion,  ^ the  i-mer  formation  energy 
and 

no=Xni  (86) 

i=l 


When  the  rotational  and  vibrational-relaxation  process  of  clusters  must  be 
taken  into  account,  an  appropriate  mechanism  must  be  assumed. 

For  small  fractions  of  clusters,  the  enthalpy,  density  and  gas 
constant  can  be  approximated  with  those  of  chemically  inert  gases,  and  the 
temperature  and  velocity  distributions  along  the  flow  axis  which  axe  obtained 
for  inert  gases  can  also  be  utilized.  Thus,  Eqs.  (83)  to  (B6 ) are  easily 
solved.  With  this  method,  dimerization  in  free  jets  was  calculated  by 
Milne  et  al  (1970)  for  A r,  Golcob  et  al  (1970)  for  Ar  and  NO,  Gordon  et  al 
(1971)  for  K,  Rb  and  Cs,  Dorfeld  and  Hudson  (1973)  for  CO2,  Beylich  and 
Robben  (1974)  for  Ar,  Seary  (1975)  for  H+Hg0*  Calo  (1975)  for  H2O,  Beylich 
(1976)  for  N£,  and  Khuth  (1977)  for  Kr,  Ar  and  Ne  (Table  4) . 


As  for  the  kinetic  process,  it  is  difficult  to  establish  a quan- 
titative scheme  without  prior  knowledge  of  the  rates  of  the  clustering 
processes.  In  the  first  stage  of  expansion,  the  cluster  kinetics  is 
primarily  governed  by  a single  cluster  species.  Once  conditions  in  the  jet 
favour  the  growth  of  this  species,  formation  of  large  clusters  will  occur 
without  a significant  rate  of  cluster  dissociation.  The  size  of  the  rate- 
controlling  cluster  may  be  a function  of  the  expansion  history.  For  fast, 
low-density  expansions,  the  size  will  be  small,  and  the  cluster  will 
consist  largely  of  dimers. 

As  for  dimer  formation  from  the  vapour,  two  processes  may  be 
considered  frcnt  kinetic  theory;  a bermolecular  mechanism  of  three-body 
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collision  with  diner  destruction  by  bimole cular  collision,  and  a stabilization 
mechanism  of  intermediate  dimer- complexes  (see  Table  4)  in  a virtual- orbiting 
state  or  of  quasibovmded  states.  It  is  well  known  that  two-body  collision 
of  hard  spheres  can  lead  to  bound  dimers  only  in  exceptional  cases.  Dissocia- 
tion of  dimers  will  occur  through  collisions  with  monomers.  Other  factors  such 
as  trimer  formation  and  dissociation  will  be  assumed  to  have  a negligible 
effect  on  the  initial  kinetics.  The  major  influence  occurs  solely  through 
dimer-monomer  interactions.  For  this  process  kinetic  equations  can  be 
developed  in  a straightforward  manner  from  collision  theory. 

A termolecular  mechanism  of  dimer  formation 


\ + *1  + \ ^**2  + *1 


(87) 

is  employed  by  Milne  et  al  (1970)  for  A r,  Golcmib  et  al  (1970)  for  A r and  NO, 
Gordon  et  al  (1971)  for  K,  Rb  and  Cs,  Dorfeld  and  Hudson  (1973)  for  C02,  and 
Calo  (1975)  for  H2O.  The  forward  rate  constant  kf  may  be  taken  as  the  normal 
termolecular  collision  rate  constant.  From  kinetic  theory, 
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with  6 - two-body  collision  diameter,  io  - approach  distance  of  the  three 
molecules  (m  &)  and  p - collision  effectiveness  written  as  0 = 30exp( - ^/kT) . 
Although  internal  degrees  of  freedom  of  the  dimer  are  not  always  equilibrated, 
for  a first-approximation  treatment  the  reverse  reaction  constant  may  be 
related  to  the  forward  reaction  rate  through  the  equilibrium  constant  Kg,  on 
the  assumption  of  microscopic  reversibility.  The  rate  of  dimer  formation  is 
then 

“ar-*f  (89) 

Golamb  et  al  (1970)  predicted  dimer  concentrations  with  0 = 1 and  io  « 6. 

The  predicted  values  were  lower  than  their  measured  values.  Milne  et  al 
(1970)  also  cal cvdated  argon  dimer  concentrations  and  looked  for  a tempera- 
ture-dependent rate  coefficient  to^  predict  their  measured  values.  They 
found  that  kf  proportional  to  either  T-^  or  T~3  could  most  nearly  reproduce 
the  experiment  al  data.  Fordon  et  al  (1971)  calculated  alkalidimer  concentra- 
tions taking  into  account  vibrational  relaxation  and  heat  of  condensation  to 
compare  their  experimental  data,  adjusting  the  collision  diameter  and  the 
extent  of  vibrational  relaxation.  They  found  that  the  collision  cross-section 
was  3.6  times  the  hard- sphere  value  and  roughly  half  or  more  of  the  heat  of 
condensation  converted  to  vibrational  excitation  of  the  dimer  molecules.  Calo 
(1975)  considered  water  dimer  formation,  using  the  collision  diameter 


( 1-665 


0.283 


(90) 


where  5o  = 2.7lX  and  T^  - 506  K.  The  calculated  results  were  in  good  agreement 
with  the  experimental  data. 
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In  the  stabilization  of  intermediate  complexes,  the  dimer  formation 
occurs  in  two  steps, 


Al+Al 


fi*+B  Ag  + B 

r 

where  Ag  is  a weakly-bound  state  which  may  or  may  not  include  metas table  and 
orbiting  pairs,  and  B is  the  monomer  or  a third  body.  Knuth  (1977)  employed 
this  model  to  yield  realistic  values  of  the  collision  effectiveness  for  dimer 
formation  of  rare  gases.  He  assumed  that  the  gas  is  in  local  thermodynamic 
equilibrium  f rcm  the  source  to  a sudden-freeze  surface  the  Ideation  df  which  is  a 
function  of  the  local  reaction- rate  coefficient.  The  relaxation  process,  that 
is,  the  mole  fraction  of  rtf  iter  is  frozen  downstream  of  the  surface.  The  location 
of  the  sudden- freeze  surface  was  given  by 


where  is  the  relaxation  time  for  the  reaction  of  dimer  formation  at 

constant  enthalpy  and  pressure,  and  Je inversely  proportional  to  the  reaction- 
rate  constant.  A quantitative  relation  between  the  terminal  dimer-mole  fraction, 
source  conditions  and  the  local  recombination  coefficient  kf  was  obtained.  From 
the  experimental  data  for  Kr,  Ar  and  Ne,  the  collision  effectiveness  was  esti- 
mated as 


(-°-6s 


and  the  terminal  mole  fraction  of  dimers  was  expressed  approximately  as 

■v>  r »./  . \7/5./,i.\2/5l5/3 


where  n0  is  the  source  nuntoer  density. 

Although  the  dimer  formation  in  free- jet  expansions  can  be  predicted 
by  the  theoretical  treatments  mentioned  above  using  the  flow  features  of  inert 
gases,  the  relevant  collision- kinetics  problems  still  remain  to  be  solved.  In 
this  regard,  Beylich  (1974)  studied  the  argon  dimerization  by  using  a Lennard- 
Jones  interaction  potential  and  a kinetic  model  to  describe  the  mass  fluxes 
between  the  vibrational  states  of  the  dimer  and  the  atom  including  the  effect 
of  the  flow  field  and  tunneling,  as  well  as  the  association  due  to  binary 
collision.  However,  further  studies  in  this  direction  are  required  for  more 
complex  substances  as  well  as  caqplex  systems  such  as  binary  mixtures  and 
heterogeneous  mixtures. 
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As  for  larger  clusters  than  dimers,  the  kinetics  axe  readily  developed 
with  bimolecular- collision  theory  between  monomer  and  i-mers,  Eq.  (2) , 


d*i 
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For  given  pre-expansion  conditions,  it  may  be  assumed  that  the  clustering 
kinetics  are  mainly  governed  by  a single  cluster  species.  Since  the  monomer 
cannot  serve  as  the  rate- controlling  species  for  a termolecular  collision,  the 
smallest  possible  rate- controlling  cluster  is  a dimer.  Once  conditions  favour 
conversion  of  dimers  to  trimers,  the  formation  of  larger  clusters  will  occur. 
Dorfeld  and  Hudson  (1973b)  developed  such  a kinetic  analysis  for  CO2  clusters. 

In  an  expansion  from  room  temperature,  most  nonpolar  gases  such  as  CO2  exhibit 
a mi nl mum  in  equilibrium- cluster  mole-fraction.  The  clustering  process  was 
assumed  to  initiate  at  the  point  where  the  equilibrium  trimer  mole-fraction 
becomes  equal  to  the  dimer  mole- fraction.  Dimer  concentrations  were  predicted 
from  the  kinetic  equation  of  dimerization  with  negligible  depletion  through 
trimer  formation.  The  observed  cluster  populations  were  explained  well  by  the 
theory. 

Beylich  and  Robben  (1974)  also  studied  the  cluster  formation  of 
larger  numbers  of  argon  atoms  (>  50) . They  measured  the  distributions  by 
means  of  a time-of-flight  mass  spectrometer,  and  compared  them  with  numerical 
calculations.  The  cluster  growth  was  adequately  described  by  kinetic  theory, 
whereas  classical  steady- state  theory  could  not  predict  the  onset  of  condensa- 
tion and  the  cluster  distribution.  Beylich  (1976)  also  compared  the  theoretical 
results  with  experiments  of  condensation  in  a nitrogen  free- jet  expansion.  The 
calculations  showed  that  classical  quasi- equilibrium  nucleation  rates  could 
yield  comparable  results  only  if  an  additional  term  in  the  free  energy  of  about 
7 KT  was  added.  Also,  the  kinetic  method  of  Pouring  (1977)  and  Kanyukhov  and 
Faizulaev  (1977),  as  well  as  the  molecular  dynamics  method  of  Zurek  and  Schieve 
(1978)  can  be  extended  to  clusterization  in  free-molecular  jets  with  sufficient 
consideration  of  the  relevant  parameters. 

Another  situation  of  practical  interest  concerning  free-molecular 
jets  is  the  problem  of  the  exhaust  of  rocket  gases  into  space.  Recently 
Wu  (1975)  examined  the  possibility  of  nucleation  and  condensation  of  certain 
species  in  the  exhaust  plume  of  an  Apollo  rocket  engine  at  high  altitudes. 
Assuming  steady  macroscopic  nucleation  theory  with  the  nonisothermal  correc- 
tion (Eq.  74)  and  using  empirical  results  for  the  plume  flow  field,  the  accu- 
mulation of  condensate  was  calculated.  The  result  showed  an  expectation  of 
water-vapour  condensation  in  the  plume  downstream  of  the  nozzle,  although  the 
nonsteady  process  of  clusterization  could  delay  the  onset  of  condensation. 

For  relevant  problems,  further  work  will  be  required  for  nonsteady,  homogeneous 
and/or  heterogeneous  nucleation  and  condensation  in  continuum  and  free-molecular 
flows. 


5.  CONCLUDING  REMARKS 

A review  has  been  given  of  recent  studies  of  homogeneous  and  hetero- 
geneous nucleation  and  condensation  and  seme  of  their  associated  steady  and 
nonstationary  flows.  For  most  experiments  in  cloud  chambers , supersonic 
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nozzles  and  shock  tubes,  homogeneous  nucleation  and  condensation  are  well- 
described  by  the  so-called  classical  theory  of  the  quasi-steady  macroscopic 
drop  model  or  its  modified  theory  with  a correction  factor.  Although  the 
factor  varies  in  the  range  from  1 to  about  10^0,  once  such  a correction  is  made 
empirically,  nucleation  and  condensation  for  most  of  these  technical  problems 
can  be  predicted  with  sufficient  reliability  by  the  gasdynamics  of  diabatic 
flows  with  phase  transitions.  Here  the  condensate  can  be  assumed  as  finely 
dispersed  with  a no-slip  velocity  from  the  gas  phase.  Nucleation  and  condensa- 
tion are  noted  primarily  for  the  effects  from  the  heat  release  of  the  latent 
heat  of  vapourization.  The  correction  factor  for  the  nucleation  theory  of  the 
macroscopic-drop  model  stems  from  a fundamental  lack  of  knowledge  of  the  volume 
and  surface  free-energy  of  gnmi i clusters  and  the  Brownian  motion  of  the  cluster. 
Thus,  knowledge  of  the  actual  values  of  a ( surface  tension),  ac  (condensation 
coefficient)  and  r (correction  factor)  remains  to  be  uncovered  before  condensation 
flow  can  be  predicted  analytically  with  confidence. 

Contrary  to  the  expansion  cloud- chamber,  the  diffusion  cloud- chamber 
produces  a steady  diffusing  stream  with  a temperature  gradient.  The  stationary 
conditions  in  the  chamber  can  be  calculated  using  one-dimensional  heat  and  mass 
diffusion  relations.  Nucleation  in  ordinary  diffusion  chambers  is  predicted 
well  using  the  classical  theory  of  the  drop  model.  By  employing  supersonic 
nozzles  or  rarefaction  waves  in  a shock  tube,  nearly- stationary  Or  nonstationary 
Uko-  or  one- dimensional  invlscid  flows  can  be  obtained  where  the  expansion  rate 
is  not  large  enough  to  produce  nonsteady  nucleation  effects.  The  observed 
static  pressure  distribution  or  history  of  the  flow  along  the  nozzle  or  tube 
axis  agrees  well  with  the  one- dimensional,  inviscid  mass,  momentum  and  energy 
flow  equations  coupled  with  a growth-rate  equation  for  the  condensate  and  a 
nucleation  rate  for  the  macroscopic- drop  model  with  a suitable  correction 
factor.  Alternatively,  a fit  of  the  numerical  calculations  to  the  observed 
results  provides  a test  for  the  nucleation  theories.  An  estimation  can  also 
be  made  of  uncertain  properties  such  as  the  correction  factor,  the  condensation 
coefficient,  and  the  macroscopic  surface-tension. 

In  general,  the  experimental  results  for  a given  substance  using 
a cloud  chamber  are  found  to  be  in  agreement  with  the  results  from  supersonic- 
nozzle  and  shock- tube  flows.  However,  the  magnitude  of  the  correction  factor 
varies  over  a very  wide  range  as  noted  above.  Thus,  when  a correction  is  used 
in  conjunction  with  the  drop  model,  the  behaviour  of  the  nucleation  process 
must  be  known  a priori  so  as  to  be  able  to  estimate  the  correction  factor.  In 
order  to  predict  such  uncertain  properties  as  well  as  to  establish  a nucleation 
theory  for  small  clusters,  two  major  approaches,  quasi- steady  and  nonsteady, 
have  been  developed.  The  first  is  a statistical-mechanical  approach  (quasi- 
steady) focusing  attention  on  the  partition  functions  of  clusters  and  the  size 
distribution  of  clusters  in  an  association-dissociation  equilibrium  state. 

The  second  is  a kinetic  approach  (nonsteady)  from  molecular  dynamics  of  cluster 
formation. 


With  the  aid  of  large  high-speed  computers,  the  evaluation  of  thermo- 
dynamic properties  of  clusters  over  a wide  range  of  size  and  temperature  becomes 
feasible  using  molecular-dynamic  or  Monte  Carlo  techniques.  Although  it  requires 
a detailed  knowledge  of  molecular  or  atomic  interactions  (potentials) , the  study 
is  in  gradual  progress  from  simple  cluster  systems  such  as  argon  vapour  to 
more  complex  systems,  e.g.,  water  vapour.  However,  it  is  unlikely  that  in  the 
near  future  sufficiently  realistic  intermolecular  potentials  will  be  developed 


for  complex- cluster  systems.  The  assumed  intermole cular  potentials  must  be 
tested  and  verified  by  comparing  the  calculated  thermal  properties  with 
experiments.  In  turn,  it  means  that  more  detailed  experiments  are  required 
for  an  understanding  of  the  features  of  small  clusters.  The  ideal  experiment 
for  such  a study  would  be  one  in  which  the  process  of  clusterization  could  be 
frozen  at  an  arbitrary  size  or  stage.  This  would  provide  detailed  information 
on  clusterization  kinetics.  One  of  the  most  promising  methods  for  this  purpose 
is  the  free-molecular-jet  technique  coupled  with  mass  spectrometry.  In  this 
system,  the  effects  of  thermal  nonequilibrium  and  relaxation  time  are  extremely 
important,  and  a detailed  treatment  of  the  kinetics  is  necessary.  The  kinetics 
of  clusterization  can  be  described  by  molecular- collision  kinetics.  Experimental 
results  from  free-molecular-beam  methods  provide  tests  for  the  kinetic  model  of 
cluster  formation.  However,  only  rather  simple  processes  of  clusterization 
such  as  dimerization  have  been  studied  to  date.  In  this  manner  a gradually- 
growing  knowledge  of  cluster  formation  will  lead  to  a better  understanding  of 
small  clusters  as  well  as  the  nucleation  process  in  free-molecule  flows  such 
as  venting-vapour  flows  of  spacecraft,  rocket  exhausts  at  high  altitude,  and 
condensed-matter  formation  in  astrophysical  problems. 

Most  nucleation  and  condensation  studies  have  been  concerned  with 
homogeneous  processes  owing  to  the  simpler  modelling  aspects  required  for  a 
theoretical  understanding.  Nucleation  of  molecules  on  ions  and  other  impurities 
is  of  considerable  importance  in  actual  and  practical  problems  such  as  the 
formation  of  clouds  and  aerosols,  metallurgical  condensation  on  substrates, 
clusterization  of  excited  species  in  combustion,  and  condensation  trails  of 
rocket  exhausts.  There  exist  many  similarities  between  homogeneous  and  hetero- 
geneous-nucleation  processes.  Current  theories  of  heterogeneous  nucleation 
are  mainly  centred  around  macroscopic  concepts,  being  essentially  extensions 
of  hamogeneous-nucleation  theories.  These  theories  can  explain  macroscopic 
features  of  the  observed  phenomena.  However,  as  in  the  homogeneous  case, 
ultimately  a theoretical  explanation  must  be  phrased  in  molecular  terms.  In 
this  regard,  it  can  be  expected  that  the  free-molecular-beam  technique  will 
provide  a promising  path  towards  understanding  the  heterogeneous  nucleation 
processes. 
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Table  1 - Cooling  Rates  Obtainable  in  Various  Flows 


Table  3 - Experimental  Studies  of  Clusterization  in  Free-Molecular  Jets 
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Table  4 - Theoretical  Studies  of  Clusterization  in  Free-Molecular  Jets 


(1)  Cluster  size:  D,  1 = 2 (dimer);  C,  i > 2 

(2)  Kinetic  model:  Ter:  A^  + A^  + A1  £ A^+A^ 
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FIG.  1 MJLTISTATE  KINETICS  OF  CLUSTERIZATION. 


FIG.  2 FREE  ENERGY  OF  CLUSTER  FORMATION  (Z£i)  AND  EQUILIBRIUM  AND  STEADY- 
STATE  NUMBER  DENSITY  OF  CLUSTER  (ni>e>  ni)  AS  A FUNCTION  OF  CLUSTER 
SIZE. 
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Vapour 


FIG.  3 HETEROGENEOUS  NUCLEATION  AND  CONDENSATION  ON  A SOLID  PARTICLE  BY 

(a)  DIRECT  IMPINGEMENT  OF  MONOMERS  ON  AN  EM3RYO  AND  (b)  DESORPTION 
OF  ADSORBED  MONOMERS  TO  AN  EMBRYO;  <j,  SURFACE  FREE  ENERGY;  a,  SURFACE 
AREA;  r,  R,  radii;  0,  CONTACT  ANGLE. 
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FIG.  5 ONSET  CONDITIONS  FOR  CONDENSATION  OF  WATER  IN  VARIOUS  FLOWS;  DIFFUSION 
CLOUD  CHAMBER  (HEIST  & REISS) , SUPERSONIC -NOZZLE  FLOW  (WEGENER  & WU) , 
SHOCK-TUBE  FLOW  (BARSCHDORFF,  KALRA,  LEE) . 
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FIG.  6 DIFFUSION  CLOUD  CHAMBER . 
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CRITICAL  SUPERSATURATION  OBTAINED  BY  A DIFFUSION  CLOUD  CHAMBER 
(KATZ  ET  AL  1975,  FIG.  5;  O-XYLENE). 
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FIG.  8 SUPERSONIC  NOZZLE  FLOW. 
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FIG.  9 NUCLEATION  AND  CONDENSATION  IN  A SUPERSONIC  NOZZLE  FLOW  (WEGENER 
ET  AL,  1972);  ETHANOL-AIR,  wq  = 0.005,  T = 105,  <r  = cr^,  a = 0.7 


FIG.  10  STEADY  FRANDTL- MEYER  FLOW. 


FIG.  11 


FIG.  12 
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SUPERCOOLING  OF  WATER  VAPOUR  AS  A FUNCTION  OF  COOLING  RATE  OBSERVED 
IN  A PRANDTL-MAYER  FLOW  (SMITH),  A NON-STEADY  LONG-DUCT  EXPANSION 
FLOW  (BARRAND  & RIEUTORD) , AND  A SHOCK-TUBE  RAREFACTION-WAVE  FLOW 
(KALRA). 
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WAVE  SYSTEM  PRODUCED  IN  A SHOCK  TUBE.  THE  SKETCHES  SHOW  THE 
PRESSURE  AND  TEMPERATURE  PROFILES,  AND  THE  WAVE  POSITION  IN 
THE  TUBE  AT  TIME  t = tj.. 
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FIG.  16  SIZE  DISTRIBUTION  OF  HETEROGENEOUS  AIT KEN  NUCLEI;  THE  DOTTED  CURVE 
IS  TYPICAL  OF  ATMOSPHERIC  AIR  AND  APPROXIMATED  BY  THE  SOLID  CURVES. 
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FIG.  18  SCHEMATIC  DIAGRAM  OF  FLOW  PROPERTIES  AMD  CLUSTERIZATION  THROUGH 
RAPID  EXPANSION  IN  A FREE  JET. 
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